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I.  INTRODUCTION 


The  objective  of  this  program  is  to  obtain  basic  information  on  the  plasma  properties  of  work¬ 
ing  fluids  for  magnetoplasmadynamic  (MPD)  power  generators.  This  information  will  be  ob¬ 
tained  through  the  use  of  diagnostic  methods  with  emphasis  on  advanced  spectroscopic  techniques. 

The  experimental  approach  of  the  program  concentrates  on  a  plasma  loop  through  which  helium 
seeded  with  cesium  circulates  through  the  test  section  at  a  variable  temperature  range  of  1000 
to  2000°K.  The  flow  rate  of  the  atmospheric  pressure  helium  is  varied  between  7  to  12  gm/sec 
to  achieve  a  Mach  number  of  0.9  through  an  area  of  approximately  1.2  cm2.  Designed  specifi¬ 
cally  for  diagnostic  investigation,  the  test  section  employs  the  basic  features  common  to  MPD 
generators.  In  preparation  for  the  operation  of  the  plasma  loop,  experimental  effort  has  been 
applied  to  the  development  of  new  spectroscopic  techniques  and  methods  of  evaluating  observed 
spectra.  Investigations  will  be  made  on  the  efficiency  of  various  nonthermal  ionization  techniques. 

The  theoretical  approach  consists  of  developing  a  theoretical  model  to  describe  nonequilibrium 
plasmas  and  supplementary  theoretical  analyses  to  aid  in  design  of  improved  MPD  power 
generators.  Instabilities  occurring  in  the  plasma  flow  will  be  investigated. 

This  report  presents  the  technical  review  of  progress  during  the  period  1  May  1963  through 
31  July  1963. 
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II.  RESUME  OF  PROGRESS 


Most  of  the  components  of  the  system  have  been  manufactured.  The  final  assembly  of  the 
test  section  will  be  completed  during  August  1963.  Based  on  the  experience  of  other  closed 
loop  experiments,  unusual  care  has  been  taken  in  the  design  of  the  system.  The  final  assem¬ 
bly  of  the  diagnostic  MPD  system  is  planned  for  the  next  reporting  period. 

Further  development  of  diagnostic  methods  was  accomplished.  A  part  of  this  work  was  the 
subject  of  a  presentation,  "A  Spectroscopic  Method  Allowing  Spatial  Resolution,  "  by  Dr.  R.  T. 
Schneider  at  the  Vlth  International  Conference  on  Ionization  Phenomena  in  Gases,  in  Paris, 
France,  8 — 13  July  1963. 

Work  on  a  theoretical  model  has  been  started.  A  portion  of  this  work  was  devoted  to  plasma 
instabilities  in  a  magnetoplasmadynamic  generator.  This  effort  was  the  subject  of  a  presenta¬ 
tion,  "Convective  Instability  of  Plasma  Flow  Across  a  Magnetic  Field,  "by  Dr.  H.  E.  Wilhelm 
at  the  Vlth  International  Conference  on  Ionization  Phenomena  in  Gases,  in  Paris,  France,  8 — 13 
July  1963. 

CLOSED  LOOP  DIAGNOSTIC  DEVICE 

The  test  section  includes  the  electrodes  and  the  windows  for  spectroscopic  measurements. 

This  section  is  made  of  tantalum  and  is  housed  in  a  stainless  steel  compartment.  Successful 
welding  of  the  tantalum  components  was  done  by  Allison.  Final  assembly  of  this  section  is 
awaiting  a  few  ceramic  parts  scheduled  to  be  shipped  in  early  August. 

An  uncooled  magnet,  for  a  20%  duty  cycle,  of  Allison  design  was  manufactured  by  a  vendor. 

The  vendor's  test  results  indicated  that  on  the  third  duty  cycle  (10  minutes  on,  50  minutes  off), 
the  temperature  rise  became  constant  and  was  115°C.  This  value  is  well  within  the  protection 
limit  given  by  the  insulation.  The  vendor's  test  results  also  showed  a  flux  density  of  13,  800 
gauss  at  60  amp  with  a  limited  power  suppiy.  The  performance  curve  shows  that  saturation 
is  not  reached  at  60  amp;  therefore,  it  appears  probable  that  a  flux  density  of  15,  000  gauss 
may  be  reached  with  the  Allison  power  supply. 

Difficulties  were  encountered  in  winding  the  tungsten  coils  for  the  heater.  However,  the  manu¬ 
facturing  procedure  has  now  been  established  and  fabrication  of  the  heater  is  progressing. 
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Status  of  other  components  is  as  follows: 

•  Received 

•  Helium  compressor  (Miehle- Dexter) 

•  50-cfm  vacuum  pump  (Welsh) 

•  Pipes,  valves,  fittings,  helium  supply  manifold 

•  Helium  and  water  flow  measuring  devices 

•  Pressure  gages  and  temperature  measuring  devices  (also  calibrated) 

•  Fabricated 

•  Cesium  supply  tank 

•  Cesium  collection  tank 

•  Two-stage  filter  (charcoal  and  absolute)  compressor  inlet 

•  Primary  condenser 

•  Concrete  pad  for  compressor  and  cooler 

•  Fabrication  in  progress 

•  Two-stage  filter  (charcoal  and  absolute)  compressor  discharge 

•  Secondary  condenser 

•  First  stage  cooler 

•  Second  stage  cooler 

•  Helium  cooler 

•  He-Cs  separator 

•  Control  panel 

•  Cs  flow  meter 

•  Compressor  aftercooler 


Final  assembly  and  snakedown  of  the  system  are  planned  for  the  next  reporting  period. 


DEVELOPMENT  AND  ADAPTATION  OF  DIAGNOSTIC  METHODS 

A  theoretical  study  was  made  to  calculate  the  absolute  line  intensities  for  cesium  as  a  seeding 
substance  in  helium.  That  portion  of  the  study  which  is  now  completed  deals  with  the  equilibrium 
values  by  thermal  excitation.  Results  are  based  on  an  extensive  parametric  study  using  the 
IBM  7090. 

In  the  experimental  part  of  the  program,  a  new  spectroscopic  technique  was  developed  which 
allows  a  spatial  resolution  of  the  plasma  properties.  This  technique  maxes  it  possible  to 


segregate  and  photograph  one  component  of  an  entirely  mixed  gas.  The  same  photograph  win 
show  the  sheaths  on  electrodes  of  the  segregated  component.  Farther,  this  method  is  a  valu¬ 
able  adjustment  aid  for  line  broadening  measurements. 

For  the  investigations  in  the  MPD-device,  this  new  method  is  expected  to  give  the  spatial  dis¬ 
tribution  of  electron  temperature  and  electron  density.  Composition  of  the  sheaths  on  the 
electrodes,  it  is  hoped,  may  be  determined  through  use  of  this  technique. 

THEORETICAL  INVESTIGATION’S 

The  theoretical  investigations  have  been  concerned  with  the  stability  of  the  plasma  flow  in  MPD 
energy  converters.  Preliminary  studies  of  the  ionization  degree  of  plasmas  which  deviate 
from  thermodynamic  equilibrium  have  been  started. 

Two  convective  processes  which  lead  to  instability  have  been  investigated.  A  one-fluid  approach 
has  shown  that  the  inertia  forces  of  the  flow  cause  a  convective  instability  with  a  growth  rate 
proportional  to  the  gradient  of  the  mean  plasma  velocity  in  the  flow  direction.  A  two-fluid 
theory  has  been  used  to  derive  a  convective  instability  with  a  growth  rate  proportional  to  the 
electron-drift  velocity.  This  latter  instability  is  conditioned  by  the  electron  and  ion  streams 
making  up  the  plasma  current. 

The  investigations  indicate  that  anomalous  or  nonequilibrium  fluctuations  are  to  be  expected 
in  the  plasma  flow,  independent  of  the  critical  Reynolds  number  above  which  the  hydromagne’i 
instability  of  the  viscous  flow  field  arises. 

The  calculations  will  be  extended  to  supersonic  flows.  Another  investigation  will  be  concerned 
with  a  possible  electrostatic  convective  instability  of  the  piasma  flow. 


Ill.  THE  ALLISON  DIAGNOSTIC  CLOSED  LOOP  MPD  DEVICE 


The  test  apparatus  utilized  in  this  program  is  a  diagnostic  device  with  features  of  a  medium 
scaled,  closed  loop,  MPD  power  generator.  To  examine  the  basic  properties  of  the  plasma, 
the  system  has  been  designed  to  permit  application  of  extensive  diagnostic  methods.  This  sec¬ 
tion  presents  pertinent  design  detail  and  operational  features  of  the  system,  heater,  test  se  tic  n, 
and  diagnostic  equipment. 

THE  SYSTEM 

A  schematic  of  the  Allison  MPD  Diagnostic  Facility  is  shown  in  Figure  1.  The  vacuum  pump, 
used  to  remove  all  of  the  air  within  the  system,  is  a  23. 5-liter/sec  (50  cfm)  Welsh  pump  'Ahich 
has  an  ultimate  vacuum  capability  of  0.01  micron.  The  piping  is  arranged  so  that  vacuum  can 
be  achieved  within  three  different  areas  of  the  loop  at  the  same  time  by  proper  valve  selection. 
The  system  will  be  evacuated  to  less  than  10  microns  at  least  twice  and  on  each  occasion  re¬ 
filled  with  purified  helium  from  the  purified  helium  source.  A  manifold  of  10  bottles  of  helium 
(two  banks  of  five  bottles  each)  installed  on  the  outside  of  the  buidling  is  separately  regulated 
so  that  one  bank  of  bottles  serves  as  the  active  bank,  the  other  as  a  reserve  bank.  The  helium 
enters  a  purification  system  consisting  of  two  NaK  bubblers  (Figure  2).  The  first  bubbler  is 
heated  electrically  to  approximately  4  80’K.  The  second  bubbler  is  allowed  to  remain  at  ap¬ 
proximately  room  temperature. 

The  two  bubblers  in  combination  have  demonstrated  the  capability  (in  other  Allison  research 
activities)  to  remove  oxygen  to  less  than  one  part  per  million.  However,  the  bubbler  system 
is  relatively  small  and  would  not  purify  the  entire  flow  of  helium  that  will  pass  through  the 
MPD  device.  Consequently,  this  bubbler  system  is  used  only  to  purify  the  gas  being  admitted 
to  fill  the  system.  A  separate  purification  system  of  a  higher  flow  capability  is  used  to  pur.fy 
the  flow  circulating  within  the  loop.  The  discharge  of  the  gas  through  the  NaK  bubbler  is  di¬ 
rected  into  various  parts  of  the  system  through  four  different  pressure  regulators.  One  pres¬ 
sure  regulator  supplies  helium  to  both  sides  of  the  diaphragm  disk  of  the  back  pressure  regulat¬ 
ing  control  valve.  If  air  were  used  for  valve  control,  there  would  be  the  danger  of  air  diifusmg 
through  the  diaphragm  into  the  system. 

A  second  regulator  is  used  to  maintain  a  positive  pressure  on  the  compressor  suction  a*  al. 
times.  This  pressure  will  vary  slightly  (in  the  order  of  a  few  millimeters  of  mercury).  The 
other  regulators  No.  3  and  4  are  used  to  supply  the  helium  necessary-  to  transport  the  cesium 
from  one  cesium  container  to  the  other. 


The  helium  compressor  is  a  Miehle-Dexter  lobe  type  (Figure  3).  The  compressor  bearings 
are  mounted  outboard  to  the  actual  helium  path,  and  there  are  carbon-face  seals  between  the 
bearing  cavity,  which  is  oil  lubricated,  and  the  compressor  cavity,  which  the  vendor  claims 
to  be  entirely  free  of  oil.  This  type  of  compressor  has  been  used  successfully  by  various 
manufacturers  of  gas  purification  systems  for  dry  boxes  and  other  purposes  where  high  purity 
gases  are  required.  This  type  of  compressor  has  a  rather  low  compression  ratio  in  order  to 
operate  in  a  nonlubricated  condition;  therefore,  the  compressor  used  requires  two  stages  to 
provide  a  1.36-atm  (gage)  output  pressure.  The  high  pressure  helium  then  is  fed  through  an 
aftercooler.  The  compressor  discharge  temperature  is  approximately  420®K.  The  aftercooler 
is  a  specially  designed  heat  exchanger  in  which  care  was  taken  to  prevent  the  admission  of 
water  from  the  cooling  side  of  the  heat  exchanger  into  the  pipe  where  the  helium  flows.  The 
stainless  steel  pipe  used  to  carry  the  helium  through  this  heat  exchanger  has  been  hydrostatically 
pressure  tested  and  ultras onically  inspected  to  rigid  requirements.  The  helium  leaves  the 
aftercooler  at  approximately  310°K  and  proceeds  through  an  activated  charcoal  filter.  The 
filter  is  possibly  not  necessary;  however,  in  the  interest  of  complete  safety,  the  filter  was 
placed  in  the  circuit  in  the  event  that  the  compressor  should  permit  slight  traces  of  oil  vapor 
to  enter  the  helium  stream.  As  the  charcoal  used  is  made  from  coconut  shell,  the  problem  of 
releasing  ammonia  found  by  other  investigators  where  animal  charcoal  was  used  should  not 
exist.  The  charcoal  filter  also  includes  an  absolute  filter  to  prevent  any  charcoal  dust  from 
being  carried  into  the  helium  stream.  This  type  of  filter  combination  has  been  used  success¬ 
fully  at  Allison  for  the  removal  of  quite  large  quantities  of  oil  vapor  from  helium  and  should 
work  very  well  in  the  removal  of  traces  of  vapor. 

Immediately  after  the  charcoal  filter  is  a  tap  for  providing  cool  helium  to  the  MPD  device 
itself.  Several  small  lines  are  fed  to  the  device  for  the  purpose  of  keeping  various  areas  in  the 
device  cool.  All  the  flow  will  be  measured.  These  small  flows  enter  the  main  helium  down¬ 
stream  from  the  MPD  device  and  do  not  alter  the  flow  through  the  device  which  is  separately 
measured.  After  passing  through  the  main  flowmeter,  the  helium  enters  the  heater  and  the 
MPD  device.  The  flowmeter  is  a  standard-gas,  ASME  sharp-edged  orifice  utilizing  a  differ¬ 
ential  pressure  gage  and  a  Heise  reference  pressure  gage.  By  the  use  of  the  aftercooler,  the 
temperature  entering  the  orifice  run  will  be  reasonably  constant,  and  by  the  use  of  the  bacK- 
pressure  regulating  valve,  the  pressure  entering  the  orifice  meter  run  will  be  relatively  con¬ 
stant.  Therefore,  only  minor  corrections  will  be  required  for  the  flow  measurement.  Down¬ 
stream  of  the  flow  measurement  device  is  a  remcie-operated  throttling  valve  which  will  be  used 
to  adjust  the  flow  through  the  heater  and  the  MPD  device.  The  adjustment  of  this  valve  will  be 
by  remote  means  at  the  control  panel  in  the  vicinity  of  the  orifice  differential  pressure  gage  so 
that  the  flow  can  be  observed  at  the  same  time  the  flow  is  being  adjusted. 
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BACK  PRESSURE  REGULATING 
COWTBOL  VALVE 


Figure  1.  Magnetoplasmadynamic  Facility  Schematic 
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Immediately  behind  the  throttling  valve  is  a  small  sampling  valve.  The  helium  can  be  circulated 
in  the  loop  at  a  low  temperature  until  samples  removed  from  the  sampling  valve  indicate  that, 
oxygen,  moisture,  hydrocarbons,  and  various  other  impurities  are  removed  to  sufficiently  low 
levels  to  permit  heater  operation.  A  common  heater  is  used  to  heat  both  the  helium  and  cesium 
to  the  required  2500°K  for  the  MPD  device.  The  heater,  however,  heats  the  cesium  and  the 
helium  in  separate  channels  utilizing  the  same  heating  element.  The  heating  is  accomplished 
with  wound  tungsten  elements  which  have  a  capacity  in  excess  of  200  kw.  The  Allison  power 
supply  available  for  use  with  this  heater  has  a  rated  capacity  of  approximately  175  kw.  This 
power  supply  is  a  remotely  operated  motor  generator  set  supplying  a  d-c  voltage  to  the  heater. 

The  helium  and  the  cesium  being  heated  separately  are  then  mixed  together  at  the  exit  of  the 
heater.  The  schematic  (Figure  1)  shows  a  space  between  the  heater  and  the  MPD  device,  but, 
actually,  these  are  closely  coupled.  The  2500°K  helium-cesium  mixture  is  then  fed  to  the 
MPD  device  which  is  explained  later  in  this  report.  The  helium  will  exit  from  the  device  at 
at  a  temperature  somewhat  less  than  the  2500°K;  however,  the  cooling  system  was  designed  to 
take  care  of  the  full  helium  and  full  cesium  flow  at  a  temperature  of  25QO°K.  The  first  cooler 
was  designed  with  a  film  coefficient  sufficiently  low  that  the  actual  pipe  or  tube  wall  would  re¬ 
main  relatively  cool  even  though  the  pipe  contained  a  helium-cesium  mixture  at  2500“K. 

The  condenser  is  water  cooled.  The  condenser  (Figure  4)  consists  of  a  tube  within  a  tube  in 
which  the  helium-cesium  mixture  flows  through  the  inner  tube,  and  water,  at  a  relatively  high 
velocity,  flows  through  the  annular  area  between  the  two  tubes.  The  first  condenser  (F.gure  5) 
is  approximately  six  meters  long.  The  exit  temperature  of  the  helium-cesium  mixture  from 
this  condenser  is  estimated  to  be  approximately  450°K.  Calculations  of  pipe  wall  temperature 
show  that  the  maximum  pipe  wall  temperature  will  be  approximately  322°K.  The  cesium  will 
start  to  condense  at  approximately  4.b  meters  Irom  entry  ot  the  heat  exchanger  condenser, 
and  most  of  the  condensation  will  be  complete  before  the  heiium-cesium  mixture  exits  from 
the  end  of  the  six-meter  run.  Pipe  wall  temperatures  at  the  end  of  the  run  are  estimated  to  be 
305°K,  which  is  still  above  the  freeze  point  of  the  cesium  that  is  rolling  along  the  bottom  of  the 
pipe.  The  cesium,  of  course,  will  not  fall  immediately  to  the  bottom  of  the  pipe,  but  will  re¬ 
main  suspended.  To  separate  the  suspended  cesium,  the  mixture  of  helium  gas  and  cesium 
liquid  will  then  be  ejected  from  the  six-meter  section  into  a  helium-cesium  separator  of  the 
centrifugal  type.  This  separator  will  remove  the  bulk  liquid  that  is  suspended  in  the  helium 
and  allow  it  to  drain  to  the  first  cesium  tank.  The  remaining  helium  will  still  contain  large 
quantities  of  cesium  vapor,  and  additional  cooling  and  condensing  is  required.  Th.s  is  accomp¬ 
lished  in  thr=e  additional  cooler-condenser  stages.  All  of  the  remaining  cooler-condenser 
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stages  are  sloped  in  such  a  manner  that  as  the  cesium  condenses  it  collects  at  the  lower  portion 
of  the  inner  pipe  of  the  double-tube  heat  exchanger  and  returns  to  the  cesium  separator  and  to 
the  cesium  tank.  The  second  stage  of  heat  exchangers  has  even  lower  pipe  wall  temperatures 
due  to  the  high  velocity  of  water  and  a  relatively  low  film  coefficient  for  the  helium-cesium 
mixture.  This  means  that  some  of  the  cesium  which  has  collected  in  the  tube  will  have  a 
tendency  to  freeze  as  it  rolls  along  the  pipe  length.  Because  the  pipes  are  sloped  toward  the 
cesium  separator,  this  cesium  will  be  constantly  rolling  in  the  direction  of  higher  temperature. 
It  is  believed  that  the  slight  amount  of  frozen  cesium  will  not  effect  a  blockage  of  the  cooler  for 
a  considerable  length  of  time. 

Provision  has  been  made  to  circulate  hot  water  through  the  helium-cesium  coolers  so  that  after 
an  experiment  Is  completed,  any  cesium  that  has  been  frozen  out  on  these  relatively  cool  pipe 
walls  can  be  melted  and  returned  to  the  separator  and  then  to  the  cesium  tank.  At  the  end  of 
approximately  12  m  of  this  type  of  helium-cesium  cooler  travel,  the  estimated  cesium  vapor 
concentration  in  the  helium  has  been  reduced  to  one  part  per  million.  At  the  completion  of  ap¬ 
proximately  18  m  of  this  type  of  cooler  travel,  the  concentration  is  less  than  five  parts  per 
billion.  There  are  approximately  2  1  m  of  helium-cesium  coolers  in  this  circuit.  At  the  end 
of  the  21-m  run,  the  helium  should  be  at  a  temperature  of  less  than  311°K.  If  for  some  reason 
the  helium-cesium  separation  has  not  been  complete,  the  helium  is  then  passed  through  a 
second  activated  charcoal  and  absolute  filter-type  combination.  Any  cesium  vapor  remaining 
in  the  helium  flow  is  then  reacted  with  the  charcoal  and  prevented  from  entering  the  helium 
compressor  where  the  cesium  might  do  damage  to  the  carbon-face  seals.  This  completes  the 
main  circuit  of  the  main  flow  of  helium  through  the  loop. 

The  purification  system,  through  which  a  part  of  the  he»ium  is  continuously  circulated,  is  con¬ 
nected  to  the  loop  at  the  main  flowmeter.  The  total  flow  in  the  compressor  s  annroxima telv 
117  liters/sec  (200  cfm).  The  MPD  device  will  require  up  to  56.4  liters/sec  of  this  total  flow. 
Approximately  2.3  5  liters /sec  (5  cfm)  will  Bow  continuously  through  the  purification  system 
regardless  of  the  flow  elsewhere  in  the  circuit;  the  remaining  35.2  liters /sec  (75  cfm)  will  flow 
through  the  back  pressure  regulating  valve.  The  flow  through  the  heater  purification  train  is 
regulated  by  a  constant  flow  regulator  and  is  measured  by  a  small  rotometer-type  flowmeter. 

The  helium  passes  through  an  economizer  heat  exchanger  (Figure  6)  where  it  ,s  preheated 
slightly  before  entering  a  gettering  trap.  The  helium  .s  admitted  to  the  copper  wool  end  of 
the  trap  and  is  heated  to  approximately  920°K.  The  heiium  then  passes  through  the  copper 
wool  gettering  trap  and  enters  the  titanium  chips  trap  where  it  .s  further  heated — to  ap- 
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proximately  1120°K;  it  is  fed  to  the  economizer  heat  exchanger  where  it  gives  up  some  of 
its  heat  to  the  incoming  gas  to  the  gettering  trap.  The  gas  then  leaves  the  economizer 
heat  exchanger  at  approximately  530°K  and  enters  a  water-cooled  heat  exchanger.  This 
heat  exchanger  is  also  designed  to  prevent  the  water  vapor  from  entering  into  the  helium 
path.  This  type  of  heat  exchanger  has  been  used  successfully  —  in  fact,  the  entire  purifi¬ 
cation  train  has  been  used  on  Allison  glove  box  purification  systems.  All  helium  piping 
in  contact  with  water  has  been  ultrasonically  inspected  ar.d  hydrostatically  pressure 
checked  to  eliminate  the  possibility  of  any  flaws  which  might  permit  water  vapor  to 
enter  the  helium  stream. 

Downstream  of  the  water  cooler  is  a  molecular  sieve  utilizing  a  Linde  molecular  sieve  material. 
This  sieve  will  remove  water  vapor  in  the  stream  which  may  have  been  introduced  with  the 
original  filling,  in  addition  to  the  water  vapor  which  might  be  generated  by  the  gettering  trap. 
The  crystalline  absorbent  sieve  material  has  a  chemical  composition  of  To.  33  ±  0.05  Na2o]  • 
£l.00  AlgOgJ  •  ^2 .48  ±  0.  03  SiC^J  •  £x  HuClJ  .  It  is  chemically  stable  to  above  600°C.  In 
addition  to  removing  water  vapor,  the  sieve  absorbs  ammonia,  should  any  appear  within  the 
system.  Because  the  sieve  and  the  charcoal  do  absorb  gases,  it  is  necessary  to  ensure  a 
sufficient  time  in  a  vacuum  for  outgassing.  The  copper  wool  and  titanium  chip  gettering  trap 
may,  under  certain  conditions,  generate  water  vapor.  If  there  is  any  hydrogen  in  the  heiium, 
the  hydrogen  will  have  a  tendency  to  reduce  the  copper  oxide  and  titanium  oxide  ;n  the  trap  to 
water  vapor. 

A  fine  stainless  steel  wire  mesh  strainer  is  located  downstream  from  the  molecular  sieve  to 
prevent  any  molecular  sieve  material  from  flaking  off  and  entering  the  heiium  stream.  The 
purified  helium  is  then  admitted  into  the  discharge  side  of  the  back  pressure  regulator  valve 
and  returned  to  the  pressure  suction  line. 

The  cesium  flow  path  begins  with  the  cesium  container  No.  2.  The  cesium  is  contained  in  a 
Mine  Safety  Appliance  shipping  container  which  includes  a  zirconium  hot  trap  and  heating  sys¬ 
tem.  The  cesium  can  be  heated  in  this  container  for  the  time  necessary  to  achieve  the  re¬ 
quired  purity.  When  cesium  is  required  in  the  system,  cesium  can  be  forced  through  the  dip 
leg  to  a  heat  exchanger  by  applying  helium  pressure  to  the  shipping  container.  The  cesium  in 
the  container  will  be  at  approximately  750°K  for  the  purification  process.  Since  the  cesium 
should  be  cooled  prior  to  entering  the  200-kw  heater  to  protect  the  heater  power  leads,  the 
cesium  will  be  cooled  to  approximately  360°K  through  a  water-type  heat  exchanger.  This  heat 
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exchanger  is  different  from  the  heat  exchanger  used  for  cooling  the  helium-cesium  mixture  in 
that  the  cesium  passes  through  a  straight  tube  around  which  is  wound  a  smailer-diameter  water 
tube.  Only  if  leaks  occur  in  both  pipes  could  the  water  ever  contact  the  cesium.  The  cesium 
flow  is  measured  with  a  heat  balance  flowmeter  which  uses  a  small  Chromoiox  electric  heating 
element  to  heat  the  cesium  through  a  small  temperature  differential.  A  wattmeter  will  be  pro¬ 
vided  to  measure  power  supplied  to  the  cesium.  A  precision  differential  temperature  measuring 
device  developed  by  Allison  will  be  used  to  measure  the  temperature  increase  of  the  cesium. 

Then  the  cesium  flow  can  be  calculated. 

Two  valves  are  provided  downstream  of  the  flowmeter,  one  a  cesium  shutoff  valve  and  the 
other  a  cesium  throttling  valve.  The  throttling  valve  is  a  remote-operated  type  actuated  from 
the  control  panel.  This  valve  will  be  air-operated  since  the  underside  of  the  diaphragm  will 
be  exposed  to  the  atmosphere  and  the  air  on  the  diaphragm  will  not  have  an  opportunity  to 
diffuse  into  the  system.  The  cesium  leaving  the  helium  flow  channel  at  the  cesium  separator 
is  directed  to  a  second  shipping  container  without  hot  trap  and  high  temperature  heating  system. 
The  low  power  heater  for  the  cesium  tank  N'o.  1  is  used  only  to  prevent  the  cesium  from  freezing. 
This  container  is  capable  of  storing  the  full  quantity  of  cesium  that  is  in  the  system.  The 
cesium  collected  in  container  No.  1  may  be  either  transferred  back  to  the  supply  container  or 
allowed  to  remain.  The  cesium  inventory  within  the  supply  container  will  be  monitored  by  an 
induction  circuit  which  remotely  locates  the  position  of  a  floating  permanent  magnet.  The  con¬ 
tinual  monitoring  of  the  cesium  level  provides  a  second  method  for  measurement  of  cesium  flow 
and  indicates  when  cesium  transfer  is  necessary. 

Cesium  in  the  supply  container  will  be  heated  to  the  gettering  temperature  and  allowed  to 
undergo  purification  for  12  hours  before  a  run  is  made.  Further  purification  will  be  made  as 
necessary  to  keep  cesium  impurities  low.  The  initial  charge  of  cesium  in  the  loop  will  be 
approximately  9  kg.  The  cesium  cooler  is  used  to  decrease  the  cesium  temperature  from 
gettering  temperature  low  enough  to  enter  the  flowmeter.  The  cooler  will  be  supplied  with 
temperate  water  to  prevent  the  cesium  from  freezing  in  the  heat  exchanger.  Valves  are  sup¬ 
plied  in  the  loop  so  that  maintenance  can  be  done  on  the  MPD  device  or  on  the  heater  without 
exposing  the  entire  loop  up  to  the  atmosphere.  The  throttling  valve  upstream  of  the  heater 
is  a  Grinnell-Saunders  diaphragm  valve  which  is  a  totally  enclosed,  mass  spectrometer  leak- 
tight  valve.  The  valve  at  the  upstream  side  of  the  activated  charcoal  and  absolute  filter  com¬ 
bination  in  the  compressor  suction  line  is  of  the  vacuum  type.  These  two  valves  can  be  closed 
to  isolate  the  heater  MPD  device  and  heat  exchangers  from  the  balance  of  the  loop.  The  vacuum 
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piping  is  arranged  so  that  this  isolated  section  can  be  evacuated  and  refilled  without  evacuating 
and  refilling  the  remainder  of  the  loop  at  the  same  time.  This  will  permit  a  savings  in  helium 
and  prevent  air  from  entering  large  area3  of  the  loop  when  work  is  required  on  the  device  or 
on  the  heater.  The  full  capacity  of  the  23.  5-liter /sec  pump  can  then  be  used  to  outgas  a  smaUer 
volume  of  the  loop. 

HIGH  TEMPERATURE  HEATER 

Figure  7  is  a  cross  section  of  the  heater.  The  tungsten  rod  in  the  center  of  the  heater  is  bored 
to  provide  a  heating  duct  for  the  cesium  flow.  At  the  end  of  the  tungsten  rod  is  one  end  o*  a 
tungsten  filament  which  is  wound  in  an  annular  area  between  the  large  tungsten  rod  and  the  outer 
shell  of  the  heater.  The  tungsten  filament  is  attached  to  an  electrode  at  the  position  where  cool 
helium  (310“K)  is  admitted,  thus  providing  cool  electrodes  for  attachment  of  the  power  leads. 
The  cesium  also  enters  this  cool  end  of  the  heater  and  is  at  a  temperature  of  about  SBOK. 


The  tungsten  filament  in  the  heater  is  supported  by  stand-off  types  of  insulators  which  are 
mounted  in  grooved  slots  on  the  large  tungsten  rod  extending  down  through  the  middle  of  the 
heater. 

The  tungsten  filament  is  designed  to  operate  at  power  levels  to  200  kw.  It  is  not  anticipated 
that  full  capacity  of  the  heater  will  be  required.  A  manufacturing  process  to  wind  the  heater 
coil  has  been  established.  Fabricated  parts  are  shown  in  Figure  8. 

TEST  SECTION 

The  MPD  test  section  has  been  designed  specifically  for  diagnostic  investigation  of  plasma 
properties.  Emphasis  has  been  placed  on  versatility,  ease  of  maintenance,  and  simplicity. 

As  shown  in  Figure  9,  the  test  section  consists  of  a  flow  channel,  electrode  blocks,  obser%'a- 
tion  ports,  and  electromagnet  assembly. 

As  the  helium-seeded  plasma  enters  the  test  section,  it  is  expanded  through  a  replaceable 
tungsten  nozzle.  By  varying  the  size  of  the  nozzle,  a  wide  variation  of  flow  rate  may  be 
achieved  for  any  given  fluid  velocity.  The  tantalum  channel  is  lined  with  high  purity  alumina 
panels  to  provide  electrical  insulation  with  respect  to  currents  induced  in  the  V  X  B  interaction. 
The  flow  stream  from  the  nozzle  is  in  the  form  of  a  free  jet  and  does  not  contact  the  channel 
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Figure  7.  MPD  Gas  Heater  Assembly 


walls  upstream  of  the  electrode  blocks.  This  feature  reduces  the  possibility  of  channel  failure 
which  could  result  from  thermal  degradation.  To  protect  the  tantalum  from  oxidation,  the  sec¬ 
tion  is  housed  within  a  1/2-in.  thick  stainless  steel  structure.  In  addition  to  providing  support 
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within  this  cavity  produces  a  slight  in-leakage  of  cool  helium  to  the  test  section. 


The  electrode  assembly  features  variable  spacing  through  the  use  of  a  bellows  sea.ing  arrange¬ 
ment.  Adjustment  of  the  two  jack-screws  provides  variation  of  the  aspect  ratio.  The  electrode 
blocks  may  be  interchanged  to  allow  investigation  of  electrode  geometry  on  plasma  properties. 

One  of  the  electrode  blocks  being  fabricated  consists  of  47  platinum  pins  shrunk  into  an  alumina 
block.  An  additional  pin  of  a  platinum- rhodium  alloy  is  inserted  into  the  alumina  block  to  allow 
thermometry  monitoring  of  the  test  section.  Electrical  leads  from  each  pin  will  be  connected 
to  a  50-pin  hermetically  sealed  connector.  External  breadboard  circu.try  permits  the  test 

7 


1 

i 


20 


vil 


I)  I) 


2/7  z  z  5  D*S  JHKU  * 

32  HQL£3  LQCATLO 
tiiTMiN  010  OTP 


Him 

11  U  Bf^aj  1 

Wll.iliMi  s 

r^Slw.  '  3$  iP||||l3ifl (I  IV  ■  ■  b^m 

■1  •>,  ,(• 

M  M 

1  pf,--  ;; : ,  .  ;f%j||§g  r«~' 

■ 

m  iii  iiiiSF 

Hiiiiiililil 

KpiB 

A  75  DtA 


r  ►  r  4  roKNl  as 

c  hon  F  “  F 


4  9GOO 


W.  J  ?LO  DtA  THAU  0  *0i  C  5 

equally  spaced  n  located 


yv/TH/N  GiO  07  t 


DC  TA/L  of  dash  + e 


T  4  £  "  Q/R  THRU  DASH  50 
1  HOI  t  i  OCA  T  ID  *V/  /  MIN  0*0  DTA 

: \m  z) 


—  /  »#•- — 


/'j#  r 


»M  »• 


AN  50/0  tO  t  0^ 

£  QU'VAL  tNTfl  rPC  JO 2 J 

I 


MS  2  4*7  i  I  OR 
IQVtvALtm  (type  301) 


I  /  ;  \ 

ft  4  -f 


t  4  40  0/4  Th/TU 


too  O'A  7 MAO 


Dt  T4H  Of  DASH  20 
2  PLACES  5 /OCA 
T HtC A NE  5  i  •  #60 


DC  TA/L  Of  OASH  2/ 

ST  OCA  T  h/CaN£  3  5  •  <3*0 


CHAM  Z?Zb  Z?/5  *45° 


f  rv  * ; 


//#  POSITION 


1  r 


ju  4)  |- 

tso  4Rt— 


r-  -  (5SC  CLOSED  POSITION 


DETAIL  Of  DASH  *4 
4  A  SIZE 


DE  TAIL  OF  OASM  /  9 

STOCK  THICKNESS'  0*0 


21 


4 


A  ^ 


*/•  ' 


••  *  *  ^  ♦ 

*  *  »/•  1 

4*  *  Lis-"* 


r/-  tail  or  dash  is 

SAME  A  5  DASH  n  EXCEPT 
as  shown 


OCT  Ail  Of  DASH  ip 
SAMI  A 5  CASH  /7£*GtPT 
AS  SHOWN 


HO  T£  5 

r  DASH  .  A»)  IN  COACEN  T  H,C  ftLbUCl  * 

SChEL  40  5  T  f  PC  3/L  -  7  AIN  I  CSS  STEEL 
LASH  2 H  FlANCtE  E/FC  Hr£ LOtN*  JlN  /SOLE 
ASA  STL  PL  hi  STA'NLt.S$  S/fEL 
1  UP E  Tuan  s  /Ac  1  Casm  *.E  a.  v 


/  DASH*/  40  *  1  70  4  M  SOS  00  I  OH  0*7  WAJ  ^  ’  *»  -  >J 
)  CASHilti  fr*  HMSCOO  241  2  S4  CL  *  025  OP  WAi  i  7  HtC 


4  4.  BOi  J  SC  At  NS  A  NUTS  MUST  SC  TrP/  302  MAN 


J  PUAChA  .IO  PHPTS  6A>*  24  A.  '.'MS,  *  4  PAP  T  HO  7mA 
PMC  A tC  AH  L  H,N  c  A" 
CHAT  ANGOLA  ?£■/*/* 


i  SAMt  AS  AN  SuOLt  5 

PuRCHHSlD  /'AN / 5  .  ASH  3:  rA r  NO  A  /OCll  5 
Lash  6  7  7 PC  4  ‘  7H/Au  L  SO  5  /EC i 

EL  C  E  ON/C  5  COPP 
MR  /  WOOD  U  d. 

9  91S-I0CS  /  ose  072  NALL  T  HIGHNESS 


9  MAH  t  EAOM  AN /OS/  32 


10  PU PC  HA  51 D  EHH7  CA  ■  Vj  H  MtO  2  M 
CANNON  Cm  EC  I  AC  CC 
32  OS  HijMIOGLljl  0/ 

LOS  A*<+CLCS  *l  AC/E 


Ft  MO* l  HUPP  6  SHAPE  L^^tE  5 
f>3 

MACH /NED  SOME  ACL  5  70  EC  v  U  0  5 


AC LD  POOS  TOOL  AM5  S7&2  AMS  3/283 
RMS  SI £34  tk.  AMS  SiQ-S 


fU  LCP  MATERIAL  TO  PE  SAME  COM  EOS/  T  /ON 
AS  AMS  7P49  (TANTALUM) 


EDM  5  T AN O AND  5/ZE  PANTS  &  MATCP/AL  (EMl/S, 

BOL  /  S.  SC  AC  *YS  ETC).  5  7ANOAAO  COM  ME  AC  /ALL  i 
AVAILABLE  FA  ATS  A  AC  ACCEPTABLE  AN  MS  OA 
ALL  /SON  S  TANDAAOS  A  AC  DC  5/PE  A  ALL  BUT  CNL  i 
/f  AFAQ/LE  AFA/LABLE 


•V  ■  + 


.  .  ■  . 

*/•  ,  V 

4f  •  J  i  -f 

4  V*  4 

*V  |  -v* 


T'/  7/7/t  OF  DASH  /S 
5 AMt  As  DASH  t  7£*C£PT 
NS  SHOW* 


5 


2*6>  OA*  THRU  N  HO/  C  s 
CQUALL  y  SP/3CCD  A  UOCA1  £  D 
WITHIN  O/O  O  T  r 


or  dash  40 


1  1 


4 — L 


M  |  W 

to  oo  q>a  thru 


DC  7 AH  O  0*5*  2  / 

5  700  Ih/ChNCSS  OtoO 


lj  '  -r 


//o  opln  position 

080  CLOSCD  POSH /ON 


or  tail  or  DASH  tof 

4  A  5/ZC 


*?•  TS"  — ■  •  J|£* 

■'  •  ; 
trrmSm 

Fl  i  ?  fl?Jmmr'2yr 


\}1‘  '  .; 

]577  Hbft  t 
*  PL  4Mf  5  5  <5  J 


—  *  t  V  • 

•  /<•<  * 


*  J,  * 

-4  — .  F-  IT  . 

SCC  Notl  2  ; 
SCI  VO  re  j  J 
'  77  NOTC  3  ] 
Ml  ts  An  io  4oi ^ 

ill  f*>ir 

t 

■\  r 

4  r  of  //Asp  t 

9  Mil  t  scrips 

« r  * 

/«  A  Mi  704  9 

I/O 

?.  < 

[7a  AMS  7 94  9  ! 

2  70 

jO  / 

T.  /S/75  1B11  J 

•  0 

*  2 

‘ r.  »m-  ye*i 

0*0 

2.  > 

‘  1  , 

0*0 

/X  / 

U  *ni  7a*  > 

CIO 

• 

i  -*>"  «■" ; 

otoo 

,  >0/V>  ♦ 

;b0 

4  , 

4-"'  "*4I  • 

C/SO 

7a  A*4S  7f*44\ 

04  O 

ff  i 

‘fvAVS  7 944 

Cbb 

’AH  l  <  i  * 

S/C  Hi  D 

01  tail  Of  OA  SH  A 

5**?/  *b  l  ASH  /7  £  •*L'*r 

AS  SHC/P/N 
AO  Tl  5 

t  DASH  *  4  -  •  //V  CONCfN  rr,C  A  /  V 

50/0  76.5  7/*/  ?//.  5TA/P/C  5  5  STCCL 
DASH  £0  fi  M'fj/.  f'/FL  HfLC/Hn  3/A  /SQlC 
ASA  b^L  TfPL  ?l*  SrA/NL  4  55  3  //'/*/ 

7tf<»4'  7  op  AS  /M  /  Of  tS  *//<■£  a/ 

/  DA’  H  •  t  AH  9  5  70  49  i  SO*  *  0*1  0*1  *•*!  -  ’  **  "»-•-' 

J  CASHtoZd  4.  4  AAJS  <f/0  274  21400*  02S0H**it  Th/C/1 

4  4i  i  £50  i  t  ■  SCALPS  3  HOTS  MUS  7  li£  7  tH£  3*2  MA7/R/A 

5  PURChA  to  PARTS  £*'■«. '*  A  "*»/,  "  AWA  *  4  A  '» 4 

AA "£/>/CAN  4.  A*  A  c  A*/* 

/  A  A  T  A  HOC  <4  A  ’4  /*// 

4  SA/**/  .75  AN  SC/OLto  5 

7  PORCH**  TO  T'AA'lO  O  4SH  J 1  '’AT  HO  A  TOC//  8 
DASH  ?h  T /PC  4  ‘ /A/A.  /  5  ’  7  f  L  i 
t<  l  A  ON/C  S  CCRP 
ma  *  wood  u  *. 

a  945  -IOC*  3  ose  07 i  HAi.L  TA/^HHLSS 

9  -PI A*,  t  f-  POP1  AH  /O  5  /  1 2 

10  PORCH  A  51 0  /HAT  CA  S  rtuOZ'P 
CANNON  Ca.CC  A  C  CO 

32  0  8  HUMBCl.  tjl  '7 
/.os  hhcllcs  7/  A/./r 


??l  pUN l  PoRR-  G  SHARP  £  -4/-  5 
6  3 

MAC  PIN  CD  SORT- ACC  5  7 0  £f  v  6*05 
*//*<?  *005  70  #7  7A?5  y;<&2  a.ms  57*7 
A  HIS  SIRA  a.  AA75  s/es 

TILL  LR  PIATT  RIAL  TO  PC  SAH)£.  COMPOSITION 
AS  ARTS  78  49  (TANTALUM) 

FOR  STANDARD  SUC  PARTS  G,  MATCPlAL  (NUTS. 
BOLTS,  SCRT  NS  £  TC  )  5  TAHOARO  COMMC  RC IALLV 
AO  AIL  ABLE  PARTS  ARC  ACCCP  TA0LC  AN  MS  OR 
ALLISON  STAROARDS  ARC  DC 5/RCAgLC  BUT  ONLY 
If  REA  OIL  Y  AfA/LABLC 


.  ZoRTC 

IL  *<•*£  * 

<  \HlA0Tf  r  ;  1 

/  >//  AO//  ^ 

^  *  '<77  A 0*7  5 

/  /  A  /  /N*Jp!  A/ AC  Ui  NOTC 

>*  hc t. ota  3/  4  NOff  4 

4  :4f$A/*i  33  1/  ;  /20 

4  Vs  »  ; 

Y  5  A//  ^S/2  HctC 


4 

2 

A  *Cf  0A*P 

- 

V' 

39  e»  Jto 

s*[ 

2 

i/7  HC7£  7 

f* 

*  j£L  NC  i  l  7 

'/•Ssr  0*  CJtlH 

i* 

* 

*  10  4 

iM 

£./>/  AM  \  132/ 

0?7 

Ys 

2 

fff  MH>\  551/' 

.120 

f ** 

40t 

4 

C pf  5  PHIS  53 1  ' 

/20 

C  S' 

/ 

ji  •  0/  LASH 

/*<', 

*3  4  4 

4t’ 

2 

Rfb  *H^  $32/ 

HU 

*09 

4  I 

7 

’(Hi  AN  .*} 

106 

4t' 

'*/  ,  4H">  J*4S 

OJ 

ft 

4HT  U  ij 

0  70 

?* 

<. 

Ah-,  12  32 

030 

t»  / 

An-.  72  72 

030 

* 

AM\  33  32 

ojo  [ 

C4 

« 

* 4  f  OP  0A1M 

<*> 

90  8  5! 

« . 

7 

U  AM '  70  4 ’/ 

060 

(*> 

J 

^  ^  AM)  70  4  9 

cto 

ft  / 

</* 

3 

[,#/  An  i>2/ 

•IQ 

P5 

V, 

! 

*  AM.,  U74 

'  -S 

•”'! 

u  \ 

* 

\a  ■  >  0*  DASH 

* 

f"'i 

\S5  59 

I 

y 

[c  Rt  AMS  53)0 

3 if  hJ’£ 

^  77/ 

ito 

7 

’/At  AM5  53/0 

c:* 

%f 

7 

PI  5  Ay*»5  9  5/0 

375 

j;;/  / 

ii 

7 

AH /0 5/32  RLf 

5/7  NO  re 

7  071 

yi’ 

7 

,  tf  AMS  >f  5 

b  ?/ 

4  , 

< 

,  /  LNC  re  /O 

I*‘ ' 

« 

'  • '  An  .  5  s  2/ 

0/0 

A3! 

*4. 

If*  ?AT5  7<97? 

QbO  \ 

>7/ 

I h  AH/OltoOto  J>77  NOT/,  t 
S  AN  SOIL  !Q  -  4  ^7£  NO*L  4 

5 /  /04toO ?  "fSLC  Note  4 
*AA  50/ 0/0  to  J>££  NC7£  4 
to  4  \ri  2  4U73  I  sec  NOT l  4 

)32  tAn ill 5 Lio  \scc  HO\c  4 

ZJJ^TJJUSAL-  ^LN£^.  4 


19 

SS'A*0 


21 


section  to  be  operated  as  either  a  segmented  or  Hall  current  generator  or  as  any  other  selected 
configuration.  With  such  a  pin-type  configuration,  the  spatial  variation  of  electrical  conduc¬ 
tivity  of  the  plasma  through  the  test  section  may  be  observed.  Further,  observance  may  be 
made  of  end  losses  which  may  develop  in  the  various  generator  types.  Since  each  electrode 
block  is  held  in  position  by  two  ceramic  setscrews,  interchangeability  is  simplified.  It  is 
anticipated  that  different  electrode  materials,  in  addition  to  various  electrode  geometries, 
will  be  investigated.  Aerodynamic  losses  may  be  determined  by  employing  blocks  with  ex¬ 
truding  pins  and  then,  for  comparison,  pins  machined  flush  with  the  electrode  block. 

Observation  windows  are  necessary  to  allow  spectroscopic  investigation  of  plasma  properties. 
Two  quartz  windows  mounted  flush  with  the  channel  and  perpendicular  to  the  electrode  blocks 
are  being  used.  The  quartz  windows  will  be  replaced  later  by  sapphire  windows.  A  continuous 
flow  of  helium  into  the  stainless  steel  structures  is  directed  over  the  windows  to  provide  cool¬ 
ing.  A  vacuum  seal  for  the  windows  consists  of  a  soft  copper  ring  seated  in  a  V-groove  on  the 
window.  To  provide  safety,  additional  quartz  windows  are  mounted  in  each  of  the  magnet  pole 
pieces. 

The  electromagnet  (Figure  10)  contains  removable  pole  pieces  which  may  be  manually  rotated 
so  that  eccentric  holes  may  be  used  to  observe  the  spatial  variation  of  plasma  properties.  The 
quartz  windows  are  mounted  within  the  core  of  the  pole  pieces  to  form  a  gas-tight  seal.  The 
gas  within  the  stainless  steel  structure  thus  provides  cooling  of  the  magnet  pole  faces.  The 
pole  pieces  protrude  through  the  stainless  steel  structure  with  an  O-ring  seal  being  used  at 
the  point  of  protrusion.  This  intermittent  duty  (10  min  on — 50  min  off)  magnet  provides  a  flux 
density  of  approximately  15,  000  gauss  with  a  gap  of  4. 45  cm. 

The  exploded  view  of  the  test  section  hardware  is  shown  in  Figure  11. 

DIAGNOSTIC  EQUIPMENT 

The  following  diagnostic  equipment  has  been  set  up,  calibrated,  and  adapted  to  the  task: 

•  A  3.4-in.  Ebert  grating  spectrograph  (Figure  12)  correcting  to  295  lines/mm  and  599  lines/ 
mm  (interchangeable)  and  having  a  theoretical  resolving  power  (first  order)  of  42,  000  and 
72,  000.  Actual  resolving  power  was  measured  to  75%  of  the  theoretical  resolving  power. 
The  instrument  was  used  to  obtain  the  spectrum  picture  camera  data  of  Section  IV  in 
this  report. 
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•  Hilger  two-prism  spectrograph  with  two  different  cameras  (Figure  13),  f/1, 5  and  f/5,  6, 
with  a  resolving  power  of  60,  000.  This  instrument  vas  used  to  detect  very  weak  cesium 
lines  for  preliminary  measurement  of  electron  densities. 

•  Microdensitometer  (Figure  14).  This  instrument  is  being  used  to  show  the  feasibility  of 
measuring  the  radial  distribution  of  broadened  cesium  lines. 

•  Spectrum  picture  camera  (Figure  15)  attached  to  a  small,  high  aperture  monochromator. 
The  spectrum  picture  camera,  described  in  Section  IV  of  this  report,  was  used  in  con¬ 
nection  with  this  small  monochromator  for  some  preliminary  studies  on  plasma  sources. 
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Figure  10.  MPD  Field  Magnet 


Figure  12.  Ebert  Grating  Spectrograph 


Figure  13.  Hilger  Two- Prism  Spectrograph 
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Figure  14.  Microdensitometer 


Figure  15.  Spectrum  Picture  Camera  Attached  to  High  Aperture  Monochromator 
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IV.  RESULTS 


In  addition  to  design  and  fabrication  of  the  diagnostic  loop,  extensive  investigations  have  been 
made  with  respect  to  development  and  adaptation  of  diagnostic  methods  and  theoretical  in¬ 
vestigations  of  plasma  instabilities  within  the  MPD  flow  channels. 


A  short  introduction  to  spectroscopic  techniques  is  presented  in  the  following  paragraphs.  This 
introduction  includes  the  basic  equations  in  this  field  and,  therefore,  provides  a  reference  for 
future  reports.  Since  the  later  reports  will  be  detailed  frequent  references  will  be  made  to  the 
relations  presented  herein. 

LINE  INTENSITY 

The  intensity  of  a  line  is  related  to  the  probability  of  transition  between  the  corresponding  en¬ 
ergy  levels. 

The  intensity  of  a  line  coming  from  an  optically  thin  plasma  is: 

I„  =/"*  i  „  d  i  U) 

*'o 

where 

i  =  radiation  power  per  unit  volume  and  unit  frequency 

ft  _  J-._4.l_  .f  -1. - 

x  -  uepui  ui  (Jictamct 

The  energy  of  a  spectral  line  is: 

i  v  *  C  Nm  (T)  h,  (2) 

where 

A™  =  transition  probability  from  level  m  to  n 
Nm  =  particle  density  of  particles  in  level  m 
h  =  Planck's  constant 
v  =  frequency 
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From  Boltzmann's  equation,  the  particle  density  is: 

N  =N0(T)i22_e'Eni/kT  (3) 

m  o  U(T) 

where 

N0  =  total  particle  density 
gm  =  statistical  weight  of  level  m 
U  =  partition  function 
Em  =  excitation  energy 
k  =  Boltzmann's  constant 
T  =  temperature 


Substitution  of  Equations  (2)  and  (3)  into  (1)  yields  the  intensity  of  the  line  per  unit  of  solid  angle: 


I  -  —  A?1  N0  (T)-^2-hv£e'Em/kT 
n  o  U(T) 


_1_ 
4  n 


(4) 


This  often-used  formula  assumes  that  the  radiation  comes  from  an  optically  thin  plasma. 


If  the  transition  probability  is  known,  Equation  (4)  can  be  used  to  determine  temperature  from 
an  absolute  measurement  of  intensity.  The  intensity  in  the  case  where  Cs  is  used  as  a  seeding 
substance  in  He  has  been  calculated  and  is  covered  in  this  report  under  subsection  Calculation 
of  the  Equilibrium  Values  of  Line  Intensities"  (also  see  Figure  19).  For  a  proper  selected 
temperature  range,  a  small  change  in  temperature  may  result  in  a  great  change  of  intensity  due 
to  the  exponential  factor.  For  those  regimes  this  method  is  very  precise. 


Equation  (4)  can  be  written: 


i  * 


In 


E 

- +  In  C 

kT 


If 


I 


=  Y, 


X, 


and  log  C  =  C1, 


(5) 


^Superscripts  denote  references  in  Section  V. 


A  linear  equation  is  obtained  with  a  slope  of  - 


1 

T‘ 


1 

Y  =  — X  +  C' 
T 


(6) 


In  using  this  method  the  intensities  of  different  lines  are  measured  and  plotted  versus  Em/k. 

If  a  straight  line  is  obtained  in  the  actual  plot,  a  criterium  is  established  for  a  Maxwell-Boltz- 

mann  distribution.  To  measure  absolute  intensity,  it  is  necessary  to  have  a  radiation  standard 

which  can  be  used  for  calibration  purposes.  In  most  cases  the  carbon  arc  with  a  temperature 

of  about  4000*K  can  be  used.  The  radiation  of  this  arc  has  been  measured  with  great  precision 
2 

by  Euler.  To  avoid  the  necessity  of  measuring  the  space  angle  and  the  absorption  in  the  glass 
optics  of  the  instrument,  it  is  advantageous  to  image  the  carbon  arc  at  the  place  of  observation 
in  the  unknown  plasma  or,  better  yet,  to  remove  the  device  containing  the  plasma  and  to  locate 
the  carbon  arc  in  the  same  place. 


EQUILIBRIUM 


Equation  (4)  is  valid  only  if  the  spectrum  line  originates  from  excited  levels  that  can  be  assumed 
to  be  populated  according  to  Boltzmann  statistics,  which  implies  a  Maxwellian  distribution  of 
electron  in  the  case  where  atoms  are  excited  mainly  by  collisions  with  the  electrons.  Further¬ 
more,  excitation  and  de-excitation  must  follow  the  same  mechanism.  De-excitation  by  spontaneous 
radiation,  therefore,  always  disturbs  Boltzmann  distribution  of  the  excited  levels  in  an  optically 
thin  plasma,  where  there  is  no  compensation  by  the  reverse  processes  of  absorption.  Thus, 
Equation  (4)  only  holds  if  de-excitation  due  to  spontaneous  radiation  can  be  neglected  compared 
to  de-excitation  due  to  collisions  of  excited  atoms  with  electrons.  According  to  Finkelnburg 
and  Maecker,  this  is  the  case  if 


+  kT  measured  in  eV)l 


where 


Ne  =  electron  number  density 

Qexc  =  collision  cross  section  for  excitation  of  level  n 
En  =  excitation  energy 

A  =  transition  probability  for  the  regarded  spectroscopic  line 


33 


LtSON 


At  electron  number  densities,  Ne  <10^,  this  condition  may  not  be  fulfilled.  Then,  the  ex¬ 
cited  levels  are  not  populated  by  Boltzmann  statistics,  and  to  relate  the  emitted  spectral 
radiation  to  plasma  parameters,  the  cross  sections  of  the  exciting  atom-electron  collisions 
must  be  known. 


RELATIVE  INTENSITIES 


If  Equation  (4)  is  applied  to  two  different  lines,  the  intensity  ratio  is 


h 

h 


Aj  g, 

A2  82 


Xh  ►  e 


-<E1  -  E2)/kT  ii 


(7) 


x  , 

The  plasma  depth  ratio  — —  is  usually  assumed  to  be  one. 

*  2 


This  is  true  only  if  the  gas  discharge 


under  consideration  has  no  temperature  profile  but  has  a  uniform  temperature.  If  a  tempera¬ 
ture  profile  exists,  the  problem  must  be  treated  with  Abel's  integral  equation  (see  subsection 
"A  Spectroscopic  Method  Allowing  Spatial  Resolution").  The  advantage  of  the  line  ratio  method 
is  that  the  particle  density  does  not  have  to  be  known. 


If  Equation  (4)  is  applied  for  the  same  line  on  a  gas  discharge  with  a  temperature  profile  on  two 
different  points: 


ix  Jx  e-(E/k)  [(1/TX)  -  (1/T0)] 

I  N„  Tj  1  (8) 

Lo  ‘N0  ux  *o 

Index  o:  point  in  core 
Index  x:  arbnrary  point 

From  this  equation  the  temperature  ratio  across  the  cross  section  of  the  discharge  can  be  estab¬ 
lished.  If  such  lines  are  selected  which  have  their  intensity  maximum  (see  Figure  19)  within 
this  temperature  range,  the  absolute  value  can  be  obtained. 

MEASUREMENT  OF  ELECTRON  DENSITY  WITH  THE  INGLIS-TELLER  TECHNIQUE 
In  the  case  of  a  hydrogen  plasma,  it  is  possible  to  obtain  a  rough  estimate  of  the  electron 

3 

density  with  the  very  simple  method  of  Inglis  and  Teller.  If  n  is  the  quantum  number  of  the 
last  visible  Balmer  line,  the  particle  density,  N,  is 

log  N  =  23.26  -  7.5  log  n  (91 
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where 


N  =  ne  +  nj  if  T  <10^/n  (ne  =  n1) 

N  =  ne  =  nj  if  T  >  lO'Vn 

Mohler  applied  this  technique  to  cesium  and  found  the  constant  to  be  23.06  rather  than  the 
aforementioned  value  of  23.26.  Using  this  method,  it  is  possible  to  determine  the  electron 
density  with  a  maximum  error  factor  of  two  only  by  observing  the  spectrum.  It  is  also  important 
to  note  that  this  method  utilizes  no  assumptions  regarding  existing  equilibrium. 

MEASUREMENT  OF  ELECTRON  DENSITY  BY  STARK  BROADENING 

The  theory  of  line  broadening  caused  by  the  Stark  Effect  has  been  extended  significantly  in  re¬ 
cent  times.  In  the  case  of  hydrogen,  the  line  profiles  given  by  Griem  and  Kolb  can  be  used. 

6  7 

For  cesium,  Agnew  and  Summes  and  Griem  have  made  calculations.  This  method  also  does 
not  depend  on  an  established  equilibrium.  It  will  be  used  quite  often  in  the  work  reported  in 
following  reports. 

CONTINUUM  RADIATION 

Free-bound  transitions  and  free-free  interactions  involve  indiscrete  energy  levels  and,  there¬ 
fore,  result  in  radiation  at  indiscrete  frequencies.  This  "continuum"  radiation  can  be  used  to 
determine  the  plasma  temperature. 


Using  Kramer's  theory,  Maecker  and  Peters  found: 


. .  .  C  CFTW  2^, 

(kT)1 


where 


GO) 


32  n  2e6 
3^3"  (2  it  n)^  I 


=  6. 36  X  10-47  cgs 


s  =  correction  factor  of  order  1. 
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DOPPLER  BROADENING 

In  the  case  of  a  low-density  plasma,  it  is  possible  to  use  Doppler  broadening  to  measure  ion 
temperature. 

/A  X 

T  =  7.  82  X  1012  M  I - 

\  * 


where 


M  =  molecular  weight 
A  X  q  =  Doppler  broadening 

In  high-density  plasmas,  the  broadening  caused  by  collision  and  by  the  Stark  effect  is  much 
larger  than  the  Doppler  broadening.  Therefore,  this  method  is,  in  general,  not  applicable  to 
high-density  plasmas. 

However,  if  an  atom  has  incompleted  shells,  such  as  the  iron  atom,  transitions  between  those 
orbits  are  almost  totally  shielded  against  Stark  interactions  and,  therefore,  allow  measurement 
of  Doppler  broadening  in  high  density  plasmas.  In  the  case  of  cesium  it  may  be  difficult  to  ob¬ 
serve  such  lines;  however,  the  ordinary  cesium  spectrum  has  some  lines — e.g. ,  the  second 
member  of  the  sharp  series  where  the  broadening  due  to  Stark  effect  is  less  than  the  Doppler 
broadening. 

DOPPLER  SHIFT 


In  the  case  of  a  high  velocity  plasma  stream,  the  measuring  of  the  Doppler  shift  of  a  spectrum 
line  would  be  a  method  to  determine  the  plasma  velocity.  The  velocity  is  given  by 


V  = 


1 


cos  a 


where 


c  =  velocity  of  light 

AX  =  Doppler  shift 

X  =  wavelength 

cos  a  =  angle  of  observation 


(12) 


^  fi 

With  a  Perot-Fabry  interferometer,  a  resolution  power - of  10  may  be  achieved.  There- 

A  X 

fore,  the  lower  limit  of  the  detectable  velocity  is  10^  cm/sec.  The  resolving  power  of  the 
Perot-Fabry  interferometer  is  given  by 


AX  X  arc  sin  G 


where 


L  =  Plate  distance 
r  -  reflectivity 


1.8  yfT 


The  useful  dispersion  is  approximately 


At  6500  A,  a  plate  distance  of  40  mm  would  be  necessary  to  get  a  resolution  power  of  about  106 
and  the  useful  spectrum  range  would  be  0.  05  A.  This  requires  a  half  width  of  the  line  to  be 
0.  01  A  or  less  for  the  line  used  in  this  type  of  investigation. 


SPATIAL  RESOLUTION 


Most  of  the  plasma  devices  under  consideration  are  cylindrical  and,  therefore,  have  a  tempera¬ 
ture  profile.  As  a  consequence,  in  order  to  apply  the  methods  previously  described,  it  is  in¬ 
sufficient  just  to  measure  the  line  intensity  or  line  broadening.  Trie  spatial  resolution  of  the 
observed  data  must  be  considered.  The  first  approach  to  get  a  spatial  resolution  is  to  form  an 
image  of  the  plasma  device  in  the  plane  of  the  entrance  slit  in  such  a  manner  that  the  axis  of 
the  cylindrical  device  is  perpendicular  to  the  entrance  slit.  Thus,  a  stigmatic  line  spectrum 
is  obtained;  the  length  of  the  spectrum  lines  gives  a  measurement  for  the  diameter  of  the  dis¬ 
charge,  and  the  intensity  distribution  along  the  line  corresponds  to  the  intensity  distribution 
across  the  diameter.  This  gives  only  the  situation  along  a  cross  section  of  infinitesimal  width. 
A  better  way  to  obtain  spatial  resolution  is  to  use  the  spectrum  picture  method,  described 
under  "A  Spectroscopic  Method  Allowing  Spatial  Resolution". 


t 


4/50*. 


TIME  RESOLUTION 

Many  plasma  devices  change  their  properties  daring  time.  These  are  pulsed  discharges  or 
fluctuating  steady-state  devices.  Both  types  are  important  in  connection  with  this  program. 
For  spectroscopic  investigations  on  such  devices,  it  is  necessary  to  employ  techniques  of 

9,10,11,12,13 

short-time  spectroscopy.  Such  techniques  are  described  in  previous  reports. 

It  is  necessary  to  have  time  resolution  as  well  as  spatial  resolution.  If  this  information  is  not 
available,  the  results  are  unambiguous  only  in  special  cases. 

Short-time  spectroscopic  techniques  will  be  employed  later  in  this  program. 
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Calculation  of  the  Equilibrium  Values  of  Line  Intensities 


To  obtain  a  more  fundamental  foundation  upon  which  to  build  the  diagnostic  experiments  and 
theoretical  model  of  the  plasma  flow  through  a  magnetic  field,  a  numerical  program  has  been 
initiated  on  the  IBM  7090  computer.  This  program  allows  the  computation  of  line  intensities 
through  which  the  observed  spectrum  data  may  be  interpreted.  Further  analyses  for  deter¬ 
mination  of  the  plasma  conductivity  have  been  provided  within  the  program  assuming  equilibrium 
at  the  electron  temperature.  The  fluid  under  investigation  is  helium  seeded  with  cesium. 

This  program,  based  on  equilibrium  conditions,  will  provide  the  reference  upon  which  to  com¬ 
pare  theoretical  and  experimental  investigations  of  nonequilbrium  plasmas.  Such  a  standard 
or  reference  is  of  extreme  importance  in  interpretation  of  observed  data. 

The  number  densities  of  ions  and  electrons  are  given  approximately  by  the  Saha  equation, 

&  .  .  .  ilW  <2-->3/2MV3/2  e-(vcs  -  A  V)/kTe  (1S> 

nCs  uCs  h3 


7 

Griem  is  of  the  opinion  that  the  error  of  number  densities  of  ions  and  electrons  given  by  the 
Saha  equation  based  on  electron  temperatures  is  10%  or  less  if 


Ne  >  9X10 


3  /kT\1/2 


17 /M  (^\ 

\eh/  \EJ 


(16) 


In  Equations  (15)  and  (16),  n^-,g+,  ne>  and  n^g  represent  the  number  densities  of  the  cesium 
ions,  electrons,  and  cesium  atoms,  respectively.  UCg  +  and  Uqs  represent  the  partition 
functions  of  the  cesium  ion  and  cesium  atoms.  E2  expresses  the  excitation  energy  of  the 
resonance  line  whereas  EH  is  the  ionization  potential  of  hydrogen.  The  ionization  potential, 
VCs’  of  the  cesium  atom  must  be  corrected  to  incorporate  the  depression,  AV,  occurring 

7 

as  a  result  of  electron  density.  The  expression  given  by  Griem  is 


r  e2 _ 

'  ~/kTe 

\  87rnee2 


=  2.  94  X  10-8, 


'V 


T„(°K) 


(17) 
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The  partition  functions  for  the  cesium  I  and  cesium  II  are  computed  on  the  basis  of  the  statistical 
weights,  gs,  and  respective  excitation  energies,  Xs,  of  the  various  energy  states. 


U, 


Cs 


gC 


-X, 


Cs 


s,  s 


s/kTe 


e'X^s  i/kT  ’^Cs  2^Te 

=  gcs,o  +  gcs,i  ’  e  +  gcs,2e  Cs'  e  +  — 

U„  +  =  V  K  +  e‘XCs.+s/kTe 

uCs  gCs,  s  e 

s  =  o 

+  +fi  +  e-^sJl/kTe  +  -XcJ2/kTe 
3,0  gCs,  le  +  gCs,2e  +-'- 


(18) 


(19) 


H 

Within  the  program,  104  energy  levels  for  cesium  I  and  63  levels  for  cesium  II  were  used  to 
compute  the  ratio  of  partition  functions  employed  in  the  Saha  equation. 


Letting  t  represent  the  right  hand  side  of  Equation  (15),  the  Saha  equation  may  be  expressed 
as: 

°Cs+  ,±  (20) 

nCs  ne 

The  total  pressure  of  the  plasma  is  given  by  the  sum  of  partial  pressures  of  the  constituents,  or 

PTOT  =  kTg  (nHe  T  °Cs'  T  kTg  (nCs+)  T  kTene  (2i) 

It  should  be  recognized  at  this  point  that  for  the  temperature  range  of  interest  in  the  proposed 
investigation  the  ionization  of  helium  is  negligible. 


By  conservation  of  charge  of  the  plasma, 
ne  =  nCs  + 


(22) 
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Letting  the  seeding  ratio,  a  ,  be  represented  as 


nCs  +  nCs+ 

a  .  - 

nHe 

and  the  ratio  of  electron  to  gas  temperature  be 

*-VTg 

the  set  of  simultaneous  Equations  (19  through  23)  may  be  solved  to  yield: 


<Ml+a+  a/J) 

/p|i(l  +  a  +  a  f)) 

2  .  PTOT  *  3 

n  +  » 

e  2(1+  a)  1 

(L  2(1+  a)  J 

T  kTe  (1  +  a) 

The  remaining  unknowns,  nCg,  nHe,  and  degree  at  ionization  are  readily  obtained. 


(23) 


(24) 


(25) 


Line  intensities  are  given  by 


where 

A™  =  transition  probability  from  level  m  to  n 
Em  =  energy  of  m-state 

Twenty-seven  of  the  strongest  cesium  lines  have  been  computed  based  on  the  transition  proba¬ 
bilities  given  in  references  15  and  16,  Should  better  values  of  transition  probabilities  become 
available,  correction  can  be  made  of  line  intensities  by  simply  multiplying  the  intensities  by 
the  ratio  of  probabilities. 

The  conductivity  of  the  plasma  has  been  computed  on  the  basis  of  calculated  electron- ion  colli¬ 
sion  cross  sections  and  assumed  values  of  electron-neutral  cross  sections.  Since  the  assumed 
values  of  cross  sections  are,  at  this  time,  questionable,  the  conductivity  values  represent 
only  order  of  magnitude  estimates.  The  calculation  has  been  specifically  programmed  so  that 
these  values  may  be  readily  changed  when  better  data  become  available. 
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LISON 


The  electron-ion  cross  section  is  given  by 

/  e2\2  [  1.5  /k3Te3\1/2 

Qi'°-  °(kTj  ln  L^(-ne  / 


(27) 


so  that  the  mean  free  path  becomes 
1 


ne  Qi  +  nCs  %  +  nHe  Qg 


(28) 


where  Qs  and  Qg  represent  the  electron-neutral  cross  sections  of  the  seed  and  gas  atoms,  re¬ 
spectively. 

The  mean  free  time  based  on  a  Maxwellian  distribution  of  velocity  is  given  as 


/8_kT\ 

\wme  / 


1/2  '  0.  62i  x  106  v/T^ 


(29) 


The  conductivity  is  then  given  as: 


<r 


0.  85  ne  e2 
me 


re  =  2.389  X  10'3  ne  re 


(30) 


Several  runs  have  been  made  with  the  program  over  a  temperature  range  of  1500  to  6000"K. 
Seeding  values  have  ranged  from  0.  5%to  100%  cesium.  The  degree  of  ionization  for  the  vari¬ 
ous  seed  ratios  is  illustrated  in  Figure  16.  As  shown  by  Figure  17,  the  optimized  seeding 
level,  based  on  assumed  cross  sections  for  the  temperature  range  of  interest,  appears  to  be 
1. 0%.  On  the  basis  of  1%  seeding  of  a  plasma  with  a  total  pressure  of  1  atm  under  conditions 
of  equilibrium  (Te  =  Tg)  the  various  properties  are  given  as  shown  in  Table  I.  For  the  same 
conditions,  the  effect  of  pressure  on  conductivity  of  pure  cesium  is  illustrated  in  Figure  18. 

As  mentioned  previously,  the  intensities  of  27  of  the  most  prominent  cesium  lines  have  been 
computed.  These  intensities  are  illustrated  graphically  in  Figure  19  and  listed  in  tabular  form 
in  Tables  II  through  VII.  With  reference  to  the  tables,  it  should  be  noted  that  only  two  signi- 

15,  16 

ficant  figures  should  be  used  since  the  values  of  transition  probabilities  are  given  to  only 
two  significant  places.  The  data  obtained  through  the  numerical  program  provides  the  basic 
information  for  interpretation  of  observed  spectra  and  further  provides  a  reference  for  the 
continued  effort  on  nonequilibrium  plasmas. 


Property 


Electron  Penalty,  cc*1 
Ion  Density  cc*1 

Cesium  Density.  cc  1 
Helium  Density,  cc*1 
Degree  of  Ionization 
Mean  Free  Path,  cm 
Mean  Free  Time.  ser 


Conductivity,  mho  m 


TABLE  I 

Theoretical  Plasma  Properties — 1%  Seeding 


Electrical  Conductivity  of  Plasmas  —  Pressure  Effect 


Figure  17.  Electrical  Conductivity  of  Plasmas 
Seeding  Effect 
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Figure  19.  Cesium  Line  Intensities 
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Cesium  Line  Intensities — (Cs  +  Cs+)/He  =0.05  (P  -  1  atm 
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Cesium  Line  Intensities — Cs  +  Cs+/He  =  0. 05  (P  -  1  atm 
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TABLE  V 


'esium  Line  Intensities — (Cs  +  Cs+)/He  =0.05  (P  -  1  atm,  Te  =  Tg) 
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esium  Line  Intensities — Cs  +  CsVHe  =  0. 05  (P  -  1  atm  -  Te  =  Tg) 
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Cesium  Line  Intensities — (Cs  +  Cs+)/He  =  0. 10  (P  =  1  atm,  1 
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43643.2617 
52692.22  17 
7652321.44 
14626994.1 
2375.66190 
16377.8989 
2C7. 568157 
1635.45288 
57.7964258 
464.  119621 

2  I.  1403039 
167.269382 


Cesium  Line  Intensities — (Cs  *  Cs*)/He  =  01.0  (P  =  1  atm,  Ti 


30003K 


3200°K 


3400’K 


3600 Tt 


3800°K 


0.539505121: 
0. 209336656 
0.1 16525746 
0. 1297 18556 
98057325.0 
0. 1081959 16 
67642329.0 
0.382453886 
0. 3876429  IE 
49234269.0 
4C8224 11.0 
1 C7C9 1  1 1 .4 
0.  162939056 
17497052.2 
0. I  10778336 
2  76  7  8999 . 7 
33441 109.7 
0.895612716 
0.178633676 
1USS  714591.773 
5635984.81 
98938.8779 
788650.523 
35769.6372 
288910.71  I 
1539 1. 7968 
122257.254 


11  0.687590706  II 
10  0.375361026  10 

12  0.151081666  12 
09  0.265167116  09 

0.200446096  09 
OV  0.204259016  09 
0.134513336  09 
09  0.760914536  09 
09  0.77C86560E  09 
0 . I C04488  76  09 
83286807.0 
22448727.2 
09  0.341639806  09 
37213760.5 
09  0.232216536  09 
60263077 . 0 
72268637.0 
09  0.156679976  10 
10  0.314171446  10 
1406677.63 
11121322.0 
2C6298. 062 
1646807.30 
77058.1514 
622848.945 
33838.7188 
2689  11.434 


0.835594  186 
0.6  16510826 
C.  186404  136 
0.488910186 
0. 369578796 
0.351081376 
G.242054C7E 
0.  1  36983486 
0.138715726 
0.  184689936 
0.  153300806 
42319159.5 
0.644  I  78846 
7 IC57202.C 
0.437762396 
0.  1  16724926 
0. I  39948906 
0.251799236 
0.507277946 
2506708.00 
19876435.7 
387081.227 
3093895.63 
148810.559 
I 2C3562 . 33 
66529.  1455 
528923.  109 


I  I  0.971207316 
10  0.936584666 
12  0.219593526 
09  0.823139616 
09  0.622230676 
09  0.555336766 
09  0.39b8l  1856 
10  0.225781636 
10  0.228550076 
09  0.31C612796 
09  0.257709266 
72666603.0 
09  0.1  1C636546 
0. 1234 12556 
09  0.751706356 
09  0.2C532076E 
09  0.246125366 
10  0.375197696 
10  0.75VC36266 
4 I0C593. 75 
32550580.0 
661911.633 
5296606.81 
261065.803 
2  112  706. 38 
1 18569.506 
943173.516 


I  1  0. I0e  I  SO  146  12 
10  0.132529016  11 
12  0.247495076  12 
09  0.  1 2  7 6 9 5 V  1 F  10 
09  0.965283656  09 
09  0.8 14734606  09 
09  0.606835376  09 
10  0.34  366  7  74F  10 
10  0.347763646  10 

09  0.481518256  09 
09  0.399249056  09 
C. 114739736  09 
10  0. 174718666  10 
09  0.196855576  09 
09  0. 118690306  10 
09  0.33  1251926  09 
09  0.397016176  09 
10  0.521 804516  10 
10  0.  105957696  1  I 
6195190.31  | 

4926  1032.5 
1041228.08 
8340377. 81 
420188.215 
3402154.94 
1936  10.984 
1 540  36 1  .  75 
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TABLE  VI 


ine  Intensities— 

■(Cs  +  Cs+)/He  =  0. 10  (P  =  1  atm,  Te  =  Tg) 

1700  K 

1800°K 

2000°K 

2200 1C 

2400 ’K 

2600  K 

2800TC 

.  1659806  IE  10 
88066.6/2 
.29057894E  10 
231V. 3451 

312. 484  62 

7195. 1736 
000.57080 
0589.5640 
1580.2393 
787.74933 
969.74451 

47. 6869b0 
1342.4783 
022.92823 
734.32452 
336.83136 
600.29234 
71622.937 
068285. 70 
07.664016 
24.3V3372 
.37224710 
7.7309561 
.78772636 
4.3116484 
.60000323 
.73771858 

0.26482560E  10 

2693196.00 

0.47627870E  10 

4  1242.5698 

31176.2346 

80363.3740 

26856.951V 

162319.543 

165372.676 

15979.8130 

13249.6025 

2603.81168 

39514.9360 

3644.32162 

26934.6279 

4885.22064 

5fc74. 99927 
1498960.70 
2624826. 88 

339. 7292 18 
2611.22318 
25.4603722 
199.030578 
6.49834222 
52.0824585 
2.25160021 

17. 7924776 

0.58095073E  10 

1 4683606.4 

0. 1 C937592E  II 
318918.949 
241078.387 
502651.277  | 

207330. 363 
1167732.39 

1  188163.50 
122932. 335 
101928.880 
21531. 08VI 
326962.316 
31323.7798 
222724. 19V 
43843.2617 
52692.2217 
7652321.44 
14626994. 1 
2375.66190 
18377.0989 
207.568157 
1635.45280 
57.7964258 

464. 1 1962 1 

2  1.  1403039 
167.269382 

0.  1 093564  7t  II 

58215024.5 

0.21 3  743  76E  II 

1682898.50 
1272142.83 
2229814.00 

1  C3C202. 25 
5808461.25 

5903856.50 
645967.852 
535601.844 
120024. 6V3 
1823607.03 
160229. 195 
1241572.31 
261338.812 
313919.594 
26749267.7 

55594136.5 

1 1545. 7404 
89709.2061 

1  143.  73390 

9040. 1 7786 
341.975155 
2750.49481 

130. 752399 

1035. 65016 

0. I6346806E  II 
0.  1816904  IE  09 
0.36996858E  II 
6665523.44 
5038627.00 

7642436.63 

3880911.63 

21900507.5 

22240628.5 
2549528.47 

2  113932.00 
497625.602 
7564051.88 
767206.937 

5  147592.  13 

1  145732.  14 
1375645.97 
05792585.0 

0. 1 675 I444E  09 

42696.6704 

333523.078 

4696.31128 

37218. 1118 

1490.05005 

12000.2372 

591.143997 

4686.39124 

0.28 142945E  II 
0.4  7  I24642E  09 
0. 58269046E  II 
21149968.0 

15987761.4 
21456445.0 

1 1803007.6 

66655543. 5 
67640372. 0 
806571 7.00 
6607659.25 

164  1267.25 
24957013.2 
2507510. 19 

16977772.5 
3961652. 16 
4754867.00 

0. 2  1 423370E  09 

0.42  1  73960E  09 

127846.328 

1002340.52 

15363.  1204 

122024.254 

5125.45508 

41324.3564 

2097.94906 

16644.2195 

0.40I6I5I0E  11 
0.IC548940E  10 

0 . 85050 120E  II 
56275  771  ,0 
42540108.0 

5  1405873.0 
30284308. 7 

0.  1 7  1 1 55  I9E  09 

0. 1 7355299E  09 

2  1406564.2 
17749173.0 
4514316.  19 

6e6665 1 2.0 
7254424.88 
46697474.0 
11346835.2 
13614393.6 

0. 4642208 IE  09 

0. 920295V5E  09 
323667. 121 
2545718.38 
41960.9912 
333920.676  ! 

14614.7789 

1  17945.489 
6144.49731 
48778.8628 

Line  Intensities- 

-(Cs  +  Cs+)/He  = 

01.0  (P  =  1  atm,  Te  =  Tg)  | 

34001C 

3600 °K 

3800°K  ' 

4000°K 

5000 ’K 

6000 1C 

8000°K 

•  8355V4  18E  11 
.6165 1082E  10 
.  1 66404  1  3E  12 
. 4889 10 1 8t  09 
. 3695 78 79fc  09 
.351081370  09 
. 242054  0 7E  09 
.136983406  10 
.  1387 15  72E  10 
. 18488993E  09 
. 1 5330080E  09 
2319159.5 
.644  1  78841:  09 
1057202.0 
.437762396  C9 
. 1 16724V2E  09 
. 1 39948 V0E  09 
.25 1 7V923E  10 
.507277946  10 
5087C8.00 
9876435.7 
87081.227 
093895.63 
48810.559 
203582.33 
6529. 1455 
28923.109 

O.97I20731E  II 
C.93658466E  10 
0.2 I959352E  12 
0.8231 398 IE  0V 
0.62223087E  09 
0. 55533676E  09 
0.39881  1 85E  09 
0.22578 163E  10 
0.22855007E  10 
0.3I08I279E  0V 

0 .257709261:  09 
72660603.0 
C.IIC63654E  10 
0.12341255E  09 
0.75  1 70635E  09 
0.2C532076E  09 
0.2  46 1 25 36E  09 
0.3  751  V769E  10 
0. 759038261:  10 
4100593.75 
32550580.0 

66191 1.633 
5296606.81 
261065.803 
2112706.38 

1  105H9.5U6 
943173.516 

0.10815014b  12 
0. 1 325290  IE  II 
0.24  7495C7E  12 
0.  12 7 695918  10 
0.96528365E  09 
0.814  73460E  09 
0.6C683537E  09 
0.343667748  10 
0.34776364E  10 
0.48  15  1 8256  09 
0 . 39924905E  09 
0.1  I473973E  09 
0.174718666  10 
0.  I9685557E  09 
0 . 1 1 869030E  10 
0.331251926  09 
0.3970I6I7E  09 
0 . 52 1 8045 1 E  10 
0. 1 C595769E  II 
6195190.31 
49261032.5 
1041228.08 
8340377.81 
420188.215 
3402154.94 
193610.984 
I54C361.75 

0.II549637E  12 
0. 1 7558055E  11 
0.267  167  76E  12 
0.  1 83  78399E  10 
0. 1 3892666E  10 
0.1 1 1 5  1  35  IE  10 

0 . 85830382E  09 
C.48622985E  10 
0.491 874  54E  10 
0.6V2  1  70  I4t  09 
0.5739  1 029E  09 

C .  1 6776062E  09 
0.255525106  10 
0.290505836  09 
0.173557426  10 
0.493865826  09 
0.591023686  09 
0.680657  I4E  10 
0. 13867943E  1  1 
8706415.38 
69334026.0 
1517451.34 

1  2  1  66 1 32. 5 
625123.687 
5063733.56 
291758. 152 
2321926.00 

0.8969955  IE  II 
0. 3094  3 I84E  II 
0.2 162378  IE  12 
0.4435772IE  10 
0. 3353 1 05  1 E  10 
0.222369776  10 
0.I9389775E  10 

0 . 1 C9970  79E  II 
0. 1  1  1 1 I842E  1 1 
f . 1 662S70CE  1 C 
0.I3787625E  10 
0.430I43I8E  09 
0.6554  760  7E  10 
0. 771 I6467E  09 
0.44495332E  10 
C .  1 3629699E  10 
0 . 1 6323669E  10 
C .  1 1  3C5545E  II 
0.23330260E  II 

19154925.5 

0. 15375I25E  09 
3040976.08 

30902353.5 
1711256.39 
13886100.0 
838574.  1  48 
6681452. 75 

0 .3690 1 36  IE  II 
0.2  198 1045E  II 
0.9  1435009E  II 
0. 3686  I042E  10 
0.29375982E  10 
0.17  15321  IE  10 
0.16254168E  10 
0.92258286E  10 
0.931 48968E  10 
0.14S23I03E  50 
0. 1 204 1 777E  10 
0.39231 182E  09 
0.59e057B5E  10 
C.7I985677E  09 
0.4058 1940E  10 
0. 1 3056958E  10 
0. 1 563 I673E  10 
0. 77203646E  10 
O.I6C680G9E  II 
15831495.9 
0.12730172E  09 
3473471.80 

280 1 0554 . 7 
1630491. 70 
13246095.5 
825387.594 
6581476.94 

0. 66386092E  10 
0.66483982E  10 

0 . 1 4459494E  II 

0. 1 5276028E  10 

0 .  1  1 54  75  1  1 E  10 
0.575I0477E  09 

0 . 6046 7  1  1 1 E  09 

0. 34354  203E  10 
0.34  652  3  37E  10 
0.56c68967c  09 
0.471 52692E  09 
0.162 1 6936E  09 
0.24  733  706E  10 
0.30632630E  09 

0.  1 6775296E  10 
0.5739  1  9  70E  09 
0.68675937E  09 
0.222272  1  7E  10 
0.46755  1  72E  10 
577  1063.63 
4663C3I 1.5 
1420629.75 

1  146944  1.0 

7  1  1848.469 
5791443.94 
3752V0.859 
2995390.50 
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A  (A) 

Ref 

for 

kaJ? 

|  Cesium  Line 

Intensities— 100%  Cesium  (P  =  latm, 

1500°K 

1600  K 

1700  K 

1800  K 

im 

I 

0«43.50 

6761.36 
6521. 10 
8079.02 
8015.71 
7609.01 
6983.49 
697'  29 
6723.28 

6586.51 
6354.98 
6217.27 

6212.87 
6034.09 
6010.33 
5844. 70 
5663.80 
4593.  18 

4555.36 

3888.37 
3876.  19 
3617.41 

3611.52 
3460.  13 

3476.88 
3400.00 
3398.00 

0.480000 

4.300000 

1.300000 
2.000000 
1.500000 
0.440000 
0.580000 

3.300000 
3.200000 
0.600000 
0.480000 
0.190000 

2.900000 
0.380000 

1.900000 
0.760000 
0.880000 
0.650000 
1.400000 
0.002900 
0.023700 
0.000938 
0.007640 
0.000550 
0.004490 
0.000320 
0.002560 

N 

R1 

BS 

JSS 

0.5876930  IE  10 

V6758I.7V7 

0.9644  708  IE  10 

7365.90210 

5568.05981 

21937.6580 

5970.00500 

33494.6167 

34212.7515 

2883.58105 

2390.91052 

406.673676 

6  163.64630 
526.806389 
4206. 75739 
647.733078 
779.973183 
614233.734 

1  125160.22 

74. 195691  1 
562.975380 
4.09092879 
31.8609645 
0.8773137  1 
7.00428343 
0.27277254 
2.14993712 

0. I0753460E  II 
3504486.66 

0. 1 8264  06  IE  11 

34672.8857 

26210.0557 

68076.470/ 

26594.831 1 

149354.230 

152408.975 

13521.2305 

1 121 I.C779 
2013.03607 
30524.8259 

2 684.46225 
20823.4636 
3409.36288 

4 103.42920 

2  1  1 7624.84 
3920577.69 

323. 877468 
2469.40045 
20.0358229 
156.496346 
4.58661 133 
36.6716800 

I. 48517881 

II. 7172027 

0. 18258 I83E  11 
10868909.0 

0.31 96420  IE  11 

135514.986 

102436.986 
299151. 746 
99007.8799 
556494. IV5 
567391.805 
52666.0947 
43667.8955 
8224.71 167 
124769.276 

1 1252.3923 
85078.9434 
14705.3827 

1  7691.5015 
6287953. V4 

1  1751332.5 

1  164.32333 
9068.1.7375 
81.0960274 
635.050789 
19.6653080 
157.432877 
6.60014212 

52.1 157470 

0.29I3I9I0E  II 

29626268.7 

0.52392624E  11 

453685.301 

342951.461 

884030.320 

317438.387 

1 785582.00 
1819167.77 

1  75784.543 
145751.100 
28643.0037 
434680.617 
40069.0435 
296292. 02  7 
53739.4448 
64627.4185 
16469186.6 
31074262.2 
3737. 16174 
28724.5334 
280.074608 

2 198.21690 
71.4644475 
572.928551 
24.7685323 
195.724604 

0. 639143706  II 

0. 16154440b  09 

0.12033194E  12 

3508645.78 

2652266.4  1 

5530010.94 

2280983.31 

12847023.7 

13071600.5 

1352462.81 

1  121389.48 

236878.256 

3597136.34 

344614.352 

245034 1.44 

482349.742 

579703.21 1 

84186429.0 

0.I6092I58E  09 

26136.2834 

202187.852 

2283.59946 

17992.7371 

635.858055 

5  106.09764 
232.578960 
1640.24504 

: 

t 

A  (A) 

.  m 

Ref 
for 
_  Am 

KAn 

Cesium  Line 

Intensities — 100  9 

Cesium  (P  =  1  at 

L. 

m, 

3000’K 

3200°K 

340CPK  i 

3600’K 

3800’K 

■ 

6943.50 
8761.38 
8521. 10 
8079.02 
8015.71 
7609.01 
6983.49 
6973.29 
6723.28 

6586.51 
6354.98 
6217.27 

6212.87 
6034.09 
6010.33 
5844.70 
5663.80 
4593.18 

4555.36 

3888.37 
3876.  19 
3617.41 

3611.52 
3480.  13 

3476.88 
3400.00 
3398.00 

0.480000 

4.300000 

1.300000 
2.000000 
1.500000 
0. 440000 
0.580000 

3.300000 
3.200000 
0.600000 
0.480000 
0.190000 
2.900000 
0.380000 

1  .900000 
0.760000 
0.880000 
0.650000 

1 .400000 
0.002900 
0.023700 
0.000938 
0.007640 
0.000550 
0.004490 
0.000320 
0.002560 

N1 

Rl 

BS 

SS 

L 

0.59995540E  12 
0.23279233E  11 
0.I29582I7E  13 
0.I4425320E  10 
0. 1090444  IE  10 
0. I203I900E  10 
0.7522I490E  09 
0.42530693E  10 
0.43 I07737E  10 
0.54 75085  IF  09 
0.45396464E  09 
0. 1 1909042E  09 
0.1611 9598E  10 
0.19457555E  09 
0.  123I9078E  10 
0.3I002775E  09 
0.371861 06E  09 
0.99596396E  10 
0. 1 98649 I4E  II 
7946600.81 
62674834.0 

1 100247.44 
8770169.25 
397775.395 
3212824.69 

171  164.105 
1359558.88 

0.77069 1 85E  12 
0.42072657E  II 
0.I6934II6E  13 
0.2972I460E  10 
0.22467 I70E  10 
0. 22894543E  10 
0. 1 507  704 1 E  10 
0.85287748E  10 
0.86403  12  IE  10 
0. 1 1 2SBB9AF  <0 
0.93352722E  09 
0.25  16 1 846E  09 
0.382929B6E  10 
0.4171 I359E  09 
0.260281 86E  10 
0.67546379E  09 
0.8 1 0029 1 3E  09 
0.I756I606E  II 
0.352I417IE  II 
15766865.1 

0. 1 2465427E  09 

2312309.28 

18458378.5 

863712.805 

6981254.94 

379284. 141 

3014116.59 

0.94756636E  12 
0.699 1 25  1 5E  II 
0.211 36285E  13 
0.55442564E  10 
0.419 1 0347E  10 
0.396 127  35E  10 
0.27449005E  10 
0. 1 553396BE  II 
0. 15730405E  II 

0 . 209665  75E  1C 
0.17  384357E  10 
0.4 7990056E  09 
0. 7305007  7E  10 
0.80579083E  09 
0.49642388E  10 

0 . 1 3236642E  10 

0 . 1 567  0248E  10 
0.28554 1 0VE  II 
0.57525475E  II 
28448825.2 
0.22539939E  09 
4389512.94 
35084871.0 

1 667516.28 
13648660. 1 
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A  Spectroscopic  Method  Allowing  Spatial  Resolution* 


The  thermionic  energy  converter  has  a  plasma  between  the  two  electrodes  which  are  only  some 
50  microns  apart.  In  addition,  there  is  an  electrode  sheath  on  both  electrodes.  Consequently, 
changes  in  the  properties  of  the  plasma  occur  within  a  very  short  distance.  It  is  highly  desirable 
to  investigate  these  changes. 


Even  if  the  plasma  is  not  confined  in  an  extremely  narrow  space,  changes  in  properties  still 
may  happen  within  short  distances.  An  example  is  the  electrode  sheaths  on  an  MPD-power 
generator. 


To  spectroscopically  obser%-e  these  changes  in  properties,  a  stigmatic  line  spectrum  is  com¬ 
monly  used.  An  image  of  the  device  under  investigation  is  formed  in  the  plane  of  the  entrance 
slit  of  the  spectrograph.  The  length  of  the  lines  obtained  corresponds  to  the  diameter  of  the 
discharge;  the  intensity  distribution  along  the  lines  is  related  to  the  intensity  distribution 
across  the  diameter.  The  true  radial  intensity-density  distribution  i(r)  can  be  found  from 
the  intensity  distribution  along  a  line,  which  may  be  called  I(x),  by  solution  of  Abel's  integral 

,  17.1B.19 

equation: 


R 

I(x)  =  2  J 

r  *  x 


i  (r)  dr 


dr 


(30) 


where  R  -  radius  of  the  discharge. 


If  a  drastic  change  in  intensity  takes  place  within  a  short  distance,  it  is  necessary  to  enlarge 
the  picture  of  the  device.  However,  enlargement  is  limited  by  the  height  of  the  entrance  slit 
of  the  spectrograph,  even  if  the  intensity  of  the  discharge  would  allow  a  greater  enlargement. 
This,  in  turn,  usually  causes  the  electrode  sheaths  to  appear  as  tiny  knots  in  the  stigmatic 
line  spectrum.  The  diffusion  effects  of  the  photographic  emulsion,  due  to  the  large  density 
gradient,  then  tend  to  falsify  the  measurement.  To  transform  the  axial  intensity  distribution 
into  radial  distribution,  the  height  or  the  width  of  the  densitometer  slit  must  be  made  ex¬ 
tremely  small  to  resolve  the  sheath  structure. 


This  paper  describes  a  method  to  avoid  these  difficulties. 


♦The  material  covered  in  this  subsection  was  presented  by  R.  T.  Schneider  as  a  paper  at  the 
VI  International  Conference  on  Ionization  Phenomena  in  Gases,  Paris,  France,  1963. 
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US  ON 


EXPERIMENTAL  SETUP 

As  stated  previously,  it  is  important  that  the  axis  of  the  cylindrical  discharge  be  parallel,  not 
perpendicular,  to  the  entrance  slit  of  the  spectrograph.  The  limitation  of  the  enlargement  by 
the  height  of  the  entrance  s lit  can  thus  be  avoided. 

The  experimental  setup  is  shown  in  Figure  20.  Axes  of  the  device  and  entrance  slit  are  parallel 
(or,  in  the  illustration,  perpendicular  to  the  paper  plane).  A  lens  forms  an  enlarged  image  of 
the  device  in  the  plane  of  the  entrance  slit.  A  slowly  rotating  mirror  is  located  in  the  light  path. 
The  rotation  of  this  mirror  causes  slow  motion  of  the  image  of  the  device  across  the  entrance 
slit.  An  exit  slit  in  the  focal  plane  of  the  spectrograph  is  located  where  a  spectrum  line  of  in¬ 
terest  appears.  Also  in  the  focal  plane  of  the  spectrograph  is  a  photographic  plate  which  can 
be  moved  slowly  in  the  direction  of  the  dispersion.  When  the  mirror  is  rotated  while  the  plate 
is  moved,  a  photographic  picture  of  the  device  is  obtained  in  the  light  of  one  specific  spectrum 
line. 


Figure  20.  Schematic  Diagram  of  Optical  Arrangement 
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A  film  instead  of  a  photographic  plate  proved  to  be  more  advantageous.  The  required  low  film 
speed  of  approximately  1  cm/minute  caused  problems  of  jerking.  Experience  showed  that  70-mm 
film  was  better  than  3  5-mm  film  in  this  respect;  therefore,  a  70-mm  camera  was  modified  to 
incorporate  a  continuous  drive  mechanism. 

The  exit  slit  is  variable,  and  the  mount  rotates  to  bring  the  exit  slit  parallel  to  the  entrance 
slit.  The  exit  slit  mount  is  attached  directly  to  the  film  window  of  the  camera  case,  reducing 
the  distance  between  the  film  and  exit  slit  to  1/10  mm. 

If  the  camera  is  used  in  connection  with  a  high  dispersion  spectrograph,  the  spectrum  lines 
are  sufficiently  separated  so  that  the  exit  slit  serves  only  as  an  aperture  and  the  image  of  the 
entrance  slit  is  used  as  an  exit  slit.  The  exit  slit  must  then  be  chosen  wider  than  the  line 
width.  The  curvature  and  the  titlt  of  the  focal  plane  can  be  neglected  over  this  short  distance. 

A  7102  photomultiplier  was  used  to  check  spectrum  line  adjustment  with  respect  to  the  exit 
slit.  A  slot  was  machined  into  the  plate  which  presses  the  film  against  the  ducts.  The  photo¬ 
multiplier  tube  was  mounted  behind  this  slot  so  that  it  looked  directly  into  the  exit  slit.  To  per¬ 
mit  this  alignment,  a  hole  is  punched  in  the  film  and  the  hole  is  then  located  in  front  of  the  slot. 
The  grating  is  rotated  until  the  spectrum  line  is  in  the  exit  slit,  as  indicated  by  obtaining  the 
maximum  signal  output  of  the  photomultiplier.  Movement  of  the  film  can  then  be  started  to 
obtain  the  spectrum  picture.  Because  the  photomultiplier  also  provides  information  concerning 
intensity,  the  proper  film  speed  can  be  selected  for  obtaining  a  correct  exposure. 

The  rotating  mirror  mechanism  was  designed  so  that  two  hours  are  required  for  one  revolution. 
Use  of  a  24-v  d-c  motor  with  a  high  gear  reduction  permits  fine  speed  adjustments  with  a  vari¬ 
able  resistor,  and  coarse  speed  adjustments  by  varying  the  focal  length  of  the  lens.  To  avoid 
distortion,  mirror  rotational  speed  is  ratioed  to  film  drive  speed  in  accordance  with  the  dif¬ 
ferent  enlargement  factors  of  the  optical  system.  This  can  be  checked  easily  with  two  marks 
placed  a  defined  distance  on  the  discharge  tube. 

APPLICATION 

This  method  was  developed  for  use  with  thermionic  energy  converters  and  MPD-power  genera¬ 
tors.  Investigations  on  both  types  will  be  reported  in  the  near  future.  This  paper,  however, 
reports  the  application  to  a  Cs  diode  which  served  as  a  checkout. 
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The  result  for  the  8521  A  cesium  line  is  shown  in  Figure  21.  The  cathode  with  an  electrode 
sheath  and  the  anode  are  visible.  Also  shown  In  Figure  21  is  the  density  distribution  of 
atoms  which  have  a  single  energy  —  for  this  case,  1.45  ev.  As  long  as  the  radiation  is  from 
an  optically  thin  layer,  the  intensity  is  given  by: 

I  =  AfJ1  h  v  n*  i  (31) 


where 


An  *  transition  probability 
h  *  Planck's  constant 
v  -  frequency 

n*  =  number  of  radiating  particles 
1  *  geometrical  sheat  thickness 


Figure  21.  Density  Distribution  of  Atoms 
Having  a  Single  Energy 
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Even  if  the  temperature  or  temperature  gradient  which  may  exist  across  the  diode  is  not  known, 
the  number  of  particles  In  the  excited  state  can  be  measured  with  this  relationship  if  (1)  the 
film  has  been  calibrated  with  the  aid  of  a  standard  lamp  and  (2)  the  transition  probability  is 
known. 

The  length  of  the  geometrical  sheath  can  be  obtained  by  using  Abel's  transformation  equation  — 
Equation  (30). 

Figure  22  is  a  comparison  of  electrode  sheaths  with  the  stigmatic  line  spectrum.  The  improve¬ 
ment  in  accuracy  can  be  readily  noted.  The  sheath  appears  in  the  line  spectrum  only  as  two 
little  knots.  In  order  to  find  the  intensity  of  the  sheaths,  the  microdensitometer  slit  has  to  be 
appropriately  small — whether  it  is  parallel  or  perpendicular  to  the  line.  In  the  case  of  Figure 
22,  however,  the  microdensitometer  slit  can  be  chosen  parallel  to  the  sheath  so  that  a  con¬ 
venient  slit  height  can  be  applied. 

For  both  the  electrode  sheath  and  stigmatic  line  methods,  the  diode  is  assumed  not  to  change 
the  intensity  with  time.  However,  if  deviations  occur,  the  error  introduced  in  the  analysis  is 
much  more  severe  for  the  stigmatic  line  spectrum. 

Figure  23  indicates  the  variations  which  accompany  a  sudden  increase  (top  of  illustration)  or 
decrease  (bottom)  in  intensity.  Using  the  picture  method,  the  variation  shows  up  as  a  bright 
or  dark  strip  across  the  diode.  The  striations  in  Figure  23  do  not  represent  a  real  geometrical 
structure.  The  device  could  have  increased  its  intensity  for  a  short  moment,  during  which  time 
the  slit  was  passing  across  the  diode  and  happened  to  be  at  the  geometrical  location  where  the 
bright  strip  now  appears.  If  .his  were  the  case,  the  stigmatic  line  spectrum  would  have  caught 
all  these  sudden  intensity  increases  which  appear  in  the  picture  as  bright  strips  and  would  have 
added  them  together.  These  strips  in  the  picture,  of  course,  make  these  parts  of  the  picture 
unusable  for  the  analysis.  However,  they  do  reveal  what  happens  to  the  device  with  respect 
to  time. 

Using  the  intensity  ratio  method,  the  electron  temperature  was  determined.  The  result  was 
3000“K  ±  20  percent.  Nine  different  intensity  ratios  were  used,  and  the  results  agreed  within 
the  experimental  error.  The  results  obtained  using  line  intensities  of  the  stigmatic  line  were 
higher  but  did  not  agree  with  each  other.  It  can  be  assumed  that  the  fluctuations  observed  in 
Figure  23  are  responsible  for  this. 


63 


L/SON 


Figure  22.  Luminous  Electrode  Sheath  as  a  Spectrum  (Upper  Part)  and  as  a 
Stigmatic  Line  Spectrum  f Lower  Part) 


*  #*■•* 
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Figure  23.  Striations  Caused  by  Short-Time  Increase  and  Decrease  of  Intensity 
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The  difficulty  of  measuring  line  profiles  is  a  disadvantage  of  the  picture  method.  But  the 
transition  from  the  picture  method  to  the  stigmatic  line  spectrum  is  easily  accomplished  by 
stopping  the  motion  of  the  rotating  mirror  and,  a  short  time  later,  the  film.  Of  course,  the 
exit  slit  has  to  be  larger  than  the  line  width.  In  Figure  22,  the  location  where  the  line  profile 
was  taken  can  be  determined  exactly.  This  profile  also  should  be  treated  with  Abel’s  transforma¬ 
tion,  Equation  (30),  as  further  described  in  Reference  20. 

If  line  radiation  of  the  nonvisible  range  is  under  consideration,  it  is  advisable  to  take  a  spectrum 
picture  before  taking  line  profiles.  This  is  a  good  way  to  determine  if  the  image  of  the  device 
is  sharp  in  the  entrance  slit  plane.  It  is  very  difficult  to  make  this  judgment  from  the  appear¬ 
ance  of  the  line  alone.  The  picture  of  the  device  in  the  plane  of  the  entrance  slit  must  be  sharp 
because  a  homogeneous  illumination  of  the  slit  would  give  only  an  average  profile. 

A  compilation  of  the  pictures  taken  by  the  described  method  is  shown  in  Figure  24.  The  pictures 
were  taken  with  the  Ebert-Grating  spectrograph.  Although  only  the  stronger  cesium  lines 
were  used,  a  considerable  change  in  the  parameters  can  be  seen.  As  the  excitation  energy 
rises,  the  electrode  sheath  density  first  increases  and  then  decreases.  This  detail  can  be  seen 
better  in  Figure  25.  The  transitions  shown  are: 

52F7/2  -  52D5/2;  72  P3 /2  -  62  Sj /2;  and  62P3/2  -  /2 

Figure  24  also  shows  that  radiation  from  higher  excitation  levels  come  more  from  the  center 
of  the  tube  due  to  the  temperature  profile  of  the  discharge.  The  tube  is  filled  with  cesium  metal. 
As  can  be  seen  in  Figure  21,  part  of  the  cesium  metal  is  evaporated;  the  rest  is  molten  and 
sticks  to  the  walls  of  the  tube.  The  anode,  illuminated  by  the  ambient  radiation,  can  be  seen 
clearly.  Because  the  cathode  is  surrounded  by  a  sheath  which  is  nearly  optically  thick,  no 
details  of  the  cathode  structure  can  be  seen. 

Figure  26  shows  the  change  of  the  cathode  sheath  due  to  variation  of  voltage  and  current  for  one 
spectrum  line.  Table  VIII  lists  the  points  on  the  iv-characteristic  (arc  mode)  where  the  pictures 
in  Figure  26  were  taken.  A  considerable  increase  in  the  sheath  thickness  with  increasing  volt¬ 
age  is  visible  in  Figure  26. 

The  described  method  has  these  features;  (1)  it  has  a  considerable  advantage  over  the  sUgmatic 
line  spectrum  if  details  of  the  discharge  are  of  interest;  (2)  it  serves  as  an  adjustment  aid  for 
line  broadening  investigations;  and  (3)  it  gives  the  exact  location  where  the  line  profile  was 


taken. 


R8 


Figure  26.  Change  in  the  Electrode  Sheath  Caused  by  Variation  of  the 
Voltage  Across  the  Tube 


TABLE  VIII 

Position  on  the  iv-Characteristic  Mode 


Point  No.  ♦ 

Volts 

Ampei 

1 

8.  1 

1 . 5 

2 

8.2 

1. 5 

3 

8.3 

1.  1 

4 

8.  5 

0.8 

5 

9.0 

0.  5 

6 

10.  0 

0.4 

♦Point  numbers  correspond  to  photo  numbers  in  Figure  26. 


In  conclusion,  the  described  method  is  advantageous  for  several  reasons: 

1.  A  single  component  of  a  mixed  gas  can  be  photographed. 

2.  In  the  case  of  nonequilibrium  investigations,  the  population  density  vs  excitation  energy 
can  be  determined  more  precisely.  This  is  true  because  a  suitable  location  can  be 
selected,  avoiding  inhomogeneities  like  sheaths  and  dark  spaces. 

3.  The  composition  of  the  sheath  can  be  determined  more  precisely. 

4.  Occasional  sudden  changes  in  intensity  can  be  detected. 

5.  The  summation  of  the  undesired  intensities  can  be  avoided. 
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Magnetohydrodynamic  Instability  of  Viscous  Plasma  Flow 

As  is  well  known  from  ordinary  hydrodynamics,  the  laminar  viscous  velocity  profile  of  a  fluid 
becomes  instable  when  in  the  perturbed  state  of  the  flow  the  inertia  force  density  is  large  com¬ 
pared  to  the  viscous  force  density.  This  condition  becomes: 


P0  v  )  v  »  MV2  v 


where  f>  =  density  field,  v  =  velocity  field,  and  m  =  viscosity,  where  quantities  referring 
to  the  unperturbed  flow  are  designated  by  zero  index  and  perturbations  by  tilde.  As  a  result 
of  this  instability,  an  anomalous  fluctuating  or  turbulent  flow  is  observed.  According  to 
Equation  (32),  the  onset  of  turbulence  can  be  expressed  in  terms  of  a  critical  Reynolds  number, 
Recr>  its  numerical  value  being  of  order  of  magnitude21 


■m. 


=  2500 


In  Equation  (33),  D  =  the  characteristic  dimension  of  the  channel;  spatial  mean  values  over  this 
length  are  marked  by  a  bar.  The  actual  range  of  Recr  varies  up  to  =s  50,  000  due  to  the  de¬ 
pendence  on  the  conditions  prevailing  in  the  channel,  in  particular,  on  the  magnitude  of  the 
initial  perturbations  (Recr  increases  strongly  for  decreasing  perturbation  amplitudes). 

in  the  case  of  an  ionized  gas,  it  is  to  be  expected  that  the  viscous  flow  instability  is  remarkably 
influenced  by  a  strong  magnetic  field.  With  regard  to  MPD  energy  converters,  plasma  flows 
with  constant  magnetic  fields  BQ  applied  perpendicular  to  the  velocity  field  vQ  (BQ  _L  vG)  are  of 
interest.  In  the  equilibrium  state  of  the  flow,  the  magnetic  field  lines  are  pulled  out  in  the  flow 
direction  in  a  proportion  corresponding  to  the  local  strength  of  the  velocity  field  — i.e.  ,  a 
magneto-viscous  boundary  layer  is  formed  near  the  channel  wall.  Its  thickness,  DB,  can  be 
calculated  approximately  by  equating  the  viscous  force  density  to  the  magnetic  induction  drag: 


y]1  [-<?0xb0>xb0' 


°B  =F 


In  Equations  (34)  and  (3  5),  the  assumption  is  included  that  the  induced  magnetic  field  is 
small  compared  to  the  exterior  magnetic  field  BQ,  which  is  valid  for  low  conductivity  plasmas 

(ix\  ~  ^<rvo  •  From  Equation  (35)  it  follows  that  Dg-*0  for  viscosity  n  -»  0  ( n 

now  denotes  the  effective  viscosity  depending  on  the  magnetic  field  in  the  plasma)”  or  electrical 
conductivity  <x  —  •=> ,  which  is  self-explanatory. 


From  Equation  (35),  it  is  recognized  that  the  effect  of  the  magnetic  field  on  the  instability  con- 
sidered  can  be  estimated  from  the  influence  of  the  magnetic  field  on  the  boundary  layer.  Re** 
placing  D  in  Equation  (33)  by  Dg  gives  the  critical  Reynolds  number,  Re®r,  for  a  plasma 
flow  perpendicular  to  a  magnetic  field: 


For  illustration,  compare  the  critical  Reynolds  numbers  for  a  plasma  flow  with  (Rec®)  and  with¬ 
out  (Recr)  magnetic  field.  From  Equations  (33)  and  (36)  follows: 


(36a) 


where  the  Hartmann  number”  M  =  D  B0\/|phas  been  introduced.  It  is  obvious  that  a  mag¬ 
netic  field  has  a  stabilizing  effect  for  large  Hartmann  numbers  Mcr  »  1  (Re®r<KRecr).  In 
this  case,  the  thickness  Dg,  of  the  magneto-viscous  boundary  layer  is  small  (gradient  of 
the  velocity  profile  only  near  the  wall)  compared  to  the  thickness  of  the  boundary  layer  of 
the  corresponding  nonmagnetoactive  plasma  flow.  For  small  magnetic  fields,  BQ-*0,  the 
presumption  upon  which  Equation  (34)  is  based  is  no  longer  valid.  Hence,  in  this  region, 
the  analysis  presented  does  not  apply. 


A  detailed  analysis  of  new  kinds  of  magnetoactive  plasma  flow  instabilities,  associated  with 
the  inhomogeneities  of  the  plasma  fields,  is  presented  in  the  following  subsection.  There, 
viscous  effects  are  neglected. 
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Convective  Instability  of  Plasma  Flow  Across  a  Magnetic  Field* 


In  analogy  to  general  convective  processes,  instabilities  associated  with  the  inhomogeneities 
of  the  plasma  are  termed  convective.  Well  known  are  the  gravitational  convective  instability  ** 
(growth  rate  «~g;  g  =  gravitational  accelera'ion)  and  the  magnetic  convective  instability” 
(growth  rate  o»~B/R;  R  =  radius  of  curvature  of  the  magnetic  field  B).  In  the  first  case,  the 
instability  is  due  to  a  transposition  of  the  perturbed  plasma  by  the  gravitational  field,  and,  in 
the  second  case,  the  instability  is  due  to  a  transposition  of  the  perturbed  plasma  by  the  tension 
of  the  magnetic  field  lines,  which  tend  to  contract  towards  its  center  of  curvature.  As  a 
general  feature,  the  growth  rates  of  convective  instabilities  are  proportional  to  the  gradients 
of  the  equilibrium  fields. 

In  various  experimental  devices— e.  g. ,  plasmadynamic  accelerators  and  energy  converters 
the  plasma  is  in  a  state  of  inhomogeneous  motion.  Under  this  condition,  an  additional  force 
density,  f  =  p  v  •  V  v,  acts  on  the  plasma  due  to  the  Inertia  of  the  mean  mass  flow.  This  force 
density  is  proportional  to  the  gradient  cf  the  mean  mass  velocity,  7,  and,  therefore,  of  the  con¬ 
vective  type.  In  the  presence  of  an  exterior  magnetic  field,  the  mean  flow  induces  an  electric 
field  which,  together  with  exterior  electric  fields,  gives  rise  to  electron  and  ion  streams  in 
the  plasma  the  difference  of  the  stream  velocities  being  proportional  to  the  curl  of  the  self- 
magnetic  field.  The  convective  instabilities  conditioned  by  the  inertia  force  of  the  mean  mass 
flow  and  by  the  induced  particle  streams  were  investigated  as  follows. 

INSTABILITY  OF  THE  MEAN  PLASMA  FLOW 

Equilibrium  State 


The  equilibrium  configuration  to  be  considered  consists  of  a  laminar,  nonviscous  plasma  flow, 
moving  with  a  nonuniform  velocity,  vQ,  paraLel  to  the  x-axis.  (x,  y,  z  are  cartesian  coordinates.) 
Perpendicular  to  the  flow  and  parallel  to  the  y-axis,  a  constant  magnetic  field,  Box,  is  applied. 

In  the  planes  z  =  t  L,  the  flow  is  bounded  by  electrodes,  which  are  connected  by  an  exterior 
circuit  containing  (1)  a  load,  R,  or  (2)  an  electric  power  supply,  P,  By  means  of  this  circuit, 
a  constant  electric  field  EQ  is  impressed  to  the  plasma. 

The  modifications  (1)  and  (2)  can,  for  example,  be  imagined  to  be  associated  with  plasma  flow 
production  by  (1)  thermal  expansion  (  'v  X  B  —  converter  ')  and  (2)  magnetic  acceleration 
("j  X  B  accelerator"). 


*The  material  covered  in  this  subsection  was  presented  as  a  paper  by  H.  E.  Wilhelm  at  the 
Vlth  International  Conference  on  Ionization  Phenomena  in  Gases,  Paris,  1963. 
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According  to  Ohm's  law,  a  current  density  is  developed  in  the  z-direction  of  the  plasma,  which 
can  be  expressed  in  terms  of  the  conductivity,  cr  ,  the  electric  field,  E*  *  Vq  X  B0,  induced 
by  the  flow  across  the  total  magnetic  field,  B0  r  B0I  +  Bq,  and  the  exterior  electric  field,  E0 
(Hall  and  electron  pressure  gradient  E?-fields  being  neglected).  The  self-magnetic  field,  Bq, 
of  this  current  is  given  by 

VXB;=%ff[Vv0X  (B01tB^)]  (37) 

with  V  •  Bq  =  0.  In  the  absence  of  viscosity,  the  velocity  field  is  a  function  of  x  alone,  = 
70(x).  From  Equation  (37)  and  V  ■  Bq  -  0,  it  follows  that  the  self-magnetic  field  is  a  function 
of  x,  Bq  =  B*  (x),  and  is  parallel  to  the  exterior  magnetic  field  Box  .  From  the  equation  of 
momentum  balance. 


p  v 
r  o  o 


7o  3  -  7Po  +  7“<V 
o 


X  B1) 


X  (B 


ox 


b*> 


it  follows  that  all  force  densities  are  parallel  to  the  flow  direction. 


(38) 


Summarizing,  the  general  solution  of  the  equilibrium  state  reads: 

P  o  =  P0  (x)>  Po  =  Po  (x)-  v0  =  v0  (x),  B*  =  bJ,  (x),  (39) 

Eq  =  constant 
tr  =  constant 
Box  =  constant 

If  the  total  electric  field  is  positive,  Eq  -  vQ  X  B0  >  0,  the  magnetic  body  force  is  negative, 

- (VXB*)XBq  <  0 — i.e.,  there  exist  solutions  with  decreasing  velocitv  field,  3  vQ/J  x<  0. 

Po 

If  the  total  electric  field  is  negative,  EQ  +  vQ  X  BQ  <  0,  the  magnetic  body  force  is  positive, 
- (  V  X  bI)  X  Bn  >0 — i.e.,  there  exist  solutions  with  increasing  velocity  field,  c)vo/dx>0. 

n.  U  U 

r0 

The  last  discussion  refers  to  the  condition  of  Figure  27 — case  1  (upper  sign),  case  2  (lower 
sign). 
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For  a  plasma  not  too  strongly  inhomogeneous  with  respect  to  the  spatial  dimension  of  the  per¬ 
turbation,  Equation  (40)  may  be  rewritten  in  the  "quasiclassical"  form 


(42) 


Another  conclusion  which  can  be  made  for  a  perturbation  of  small  spatial  extension  (  X  =  2*  /|ki-0) 
is  that  in  Equation  (40)  the  amplitude  q |k|-»0  of  the  elementary  long  wavelength  (X-*<°)  per¬ 
turbations  is  small  compared  to  the  amplitude  q[k|  —  »  of  the  elementary  short  wavelength  (  X  -0) 
perturbations.  For  this  reason,  the  following  inequality  can  be  derived  between  the  spatial 
derivatives  "v"  of  equilibrium  fields,  qQ  (?),  and  perturbations,  qn(r,  t): 

ojSn^v]  ~  t0[%  <*%>]'  ‘  <<=  1  <43> 

where  0  [x]  =  order  of  magnitude  of  x. 

For  study  of  a  macroscopic  instability,  the  macroscopic  approach  is  used.  Dissipation  and 
space  charges  are  neglected  throughout.”  The  field  equations  are  the  linearized  equation  of 
momentum  conservation: 


-  1  ~ 

V  B  +  —  B  •  V  Bq 
^o 


— • 

B  ■  B\  1 

=  -  v(p+^r)+  ^B°‘ 

and  the  linearized  equations  of  conservation  of  mass,  energy,  and  magnetic  flux  density: 

[^  +  v°  •  V  j  p  =  -  (v  •  VPo  +  7  Po  V  •  V) 

'1  +  7Q  •  V  j  B  =  -  (v  •  V  B0  +  B0  V  •  7  -  B0  •  V  7) 

In  deriving  Equations  (45),  (46),  and  (47),  the  unimportant  terms-*, 


(44) 


(4  5) 

(46) 

(47) 


\pv  ■  vG  |«  |p0  v  •  v|,  |  yP  v  •  v0j«)yp0  v  •  v|. 

|b  V  •  vQ  -  B.  V  V0|«|b0V  .  v  -  B0  •  V  v 


(4  8) 
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have  been  suppressed.  Elimination  of  P,  p,  and  B  from  Equation  (44)  by  means  of  Equations 
(45),  (46),  and  (47)  gives  the  characteristic  equation  for  the  velocity  perturbation  v: 


'a  -  l2  ~ 

‘  a  -  1 

_  +  Vo .  VJ  v  +  p0 

v°‘rJ 

V  p0  +  Pq  V  •  v)  v0  •  v  vc 


=  V  (v  ■  C  p0  +  yp0  v  '  v)  +  (v  •  V  B0  *  bq  v  •  v  -  B0  •  V  v)  •  —  - 

Po. 


—  B0  •  V  <-B0  -r  V  +  B0r  •  7)  - 
Mo 

—  (v  •  V  B0  +  B0  V  •  v)  .  C  B0 


(49) 


The  last  magnetic  convective  expression  vanishes  as  the  B0-lines  are  straight-lined. 

The  general  solution  of  Equation  (49)  is  rather  complicated.  It  is  suitable  to  specify  instabilities 
propagating  in  the  three  distinguished  directions  of  the  system. 


Perturbations  Propagating  Parallel  to  the  Flow 


For  k  =  kx,  it  follows  from  the  components  of  Equation  (49)  by  consideration  of  the  equilibrium 
condition,  Equation  (38):  vv  =  0,  vz  =  0,  and 


—  +  vo  Vx 


at 


yp0 


' 2  Q  av0  r  a  ] 

J  +  p°  TT  b +  Vo  rxJ 


^x2  vx. 


dv0  dp0) 

°  ax  ax )  Vx 


(50) 


From  Equation  (50),  the  dispersion  equation  is  derived  as  condition  for  a  nontrivial  solution: 


o<  -t-  i  kx 


1 

2 


t  4 


ain  Pp  _  VZ  Cs2  +  Cb2 

ain  vn  (Svn/dx)2 


(51) 


[o  o  1 1  /2  2  22 

Cs+Cgl  ,  where  Cs  =  y  p0/P0  and  Cg  =  B0/  P0  P0,  is  the  wave  velocity  for  magneto¬ 
acoustic  perturbations  propagating  perpendicular  to  Bq. 
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The  conditions  for  instability  (Re  o>  >0)  are  obtained  as: 


I. 


dv0 

dx 


<  0 


(52) 


II. 


dv 


din  Pn 


—  >  0  and 
ox  din  vr 


>0. 


(53) 


Against  the  perturbations  regarded,  according  to  Equation  (52),  a  decreasing  flow  field  is 
always  unstable;  according  to  Equation  (53),  an  increasing  flow  field  is  unstable  only  when  the 
condition  din  PQI  din  vQ  >  0  is  satisfied.  In  the  last  case  only  long  wavelength  instabilities, 

9  din  A  ( dvQ/d  x)2 

for  which  <  - -  .  - ,  are  excited. 

x  din  vQ  c2  +  C  2 
s  B 

Perturbations  Propagating  Parallel  to  the  Magnetic  Field 


For  k  -  ky,  it  follows  from  the  components  of  Equation  (49)  by  consideration  of  the  equilibrium 
condition,  Equation  (38): 

a2  2  b-h 


[  d2  2  B2', 

LP°  a?  -  ry  ■£]' 


0*  and 
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-  P  n  V, 
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°V°d^  Vy  ^ 


(54) 


(55) 


The  determinant  of  Equations  (54)  and  (55)  set  to  zero  leads  to  the  dispersion  equation: 


&)2  +  k2  C2 
Ky 


=  + 


(  o  dvQ  2  din  v  dlnPQ  2  ^  2\ 


(56) 


♦Solution:  o>  =  ±  iky  Cg  (Alfven-wave  propagating  parallel  to  BQ). 
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The  convective  solutions  of  Equation  (56)  read,  in  the  limiting  case  of: 


•  Short  wavelengths 


lim  <*>  =  — 

|kyb«  2 

Long  wavelengths 


i  k  CB1 
l  J  k„  =  » 


lim  01 


.  JL  ta.  j 

2  dx  ( 


1  ±1/1  +  4 


din  PQ 
din  v0 


By  applying  the  Routh-Hurwitz  criterion  to  the  polynomial  of  Equation  (56),  the  same  instability 
conditions  are  obtained  as  for  perturbations  propagating  parallel  to  the  flow. 

Perturbations  Propagating  Parallel  to  the  Current 

Forlc  ="kz,  it  follows  from  the  components  of  Equation  (49)  by  consideration  of  the  equilibrium 
condition,  Equation  (38):vy  -  0  and 


d2  dvQ  3  ^vodpol~  „  dvo  _  ~ 

0  °  3x  Jt  °  dx  dx  dx 

d2  2  /  b2\1  .  dv0  _ 

L  o  "dt^  Vzko+-)Jvz=-PoV0  ,x  -zvx 


The  determinant  of  Equations  (59)  and  (60)  set  to  zero  leads  to  the  dispersion  equation: 


’  2  +  k2  (C2 


l  +C|)  ("2  +  vo 

2  /  2  ^ln  vo\2 

k*  •  k  — ) 


din  v  0  dlnPo' 
dx  dx  i 


The  convective  solutions  of  Equation  (61)  read  in  the  limiting  case  of: 


•  Short  wavelengths 


lim  oi 

lkzl—  00 


1  avo  j 

2  dx  } 


1  +4 


din  P0  +  Vq 
din  v0  Cg  +  Cg 


•  Long  wavelengths 


By  applying  the  Routh-Hurwitz  criterion  to  the  polynomial  of  Equation  (61),  the  same  instability 
renditions  are  obtained  as  for  perturbations  propagating  parallel  to  the  flow. 

The  convective  instability  considered  is  conditioned  by  the  inertia  forces  of  the  flow  and  vanishes 
when  the  latter  becomes  uniform.  It  is  due  to  fluctuations  of  the  pressure,  density,  magnetic, 
and  velocity  fields.  Remarkable  is  the  fact  that  a  flow  slowing  down  always  becomes  convectively 
unstable. 

CONVECTIVE  INSTABILITY  OF  THE  PARTICLE  STREAMS 


Equilibrium  State 


In  this  section  special  regard  is  given  to  a  convective  process  connected  with  the  electron  and 
ion  streams,  which  could  not  be  comprehended  by  the  previously  discussed  one-fluid  approach. 
In  order  to  study  the  effect  in  pure  form,  the  inertia  forces  of  the  mean  mass  flow  are  elimin¬ 
ated  by  appropriate  modification  of  the  plasma  model  (Figure  27).  Suppose  the  channel  is 
bounded  by  additional  nonconducting  walls  perpendicular  to  the  flow,  and  that  the  current  is 
produced  by  means  of  the  exterior  power  supply.  Then,  in  the  equilibrium  state,  the  mean 
mass  velocity,  v0,  is  zero,  and  the  magnetic  body  force  is  balanced  by  the  pressure  gradient 
alone: 

0---Vpo  +  -(VX?o)X(Boi  +  Bq)  (64) 

^o 

The  difference  between  the  ion  and  electron  drift  velocities  in  the  quasineutral  plasma, 


£  nos  es  =  0  (65) 

s=e,  i 


(e,  i  index  for  the  electron  and  ion  component  respectively)  is  given  from  Maxwell's  equation 
by: 


V  X  B 


I 


(66) 
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(67) 


The  drift  velocities  are  constant.  The  total  momentum  of  the  streams  vanishes: 

E  nr 
s=e,  i 


■os  ms  v  iS  =  ^ 


Finally,  the  condition  that  all  forces  acting  on  the  plasma  components  s  =  e,  i  are  in  equilibrium 
reads: 


fW 


*ox 


0  =  *  ^  Pos  *  nos  es  |E0  *  v  xs  *  ©c 

^  nos  nor  ^  sr  is  "  ^  ir^ 
r=e,  i 


(68) 


where  fsr  is  the  friction  coefficient  of  the  s-component  (fei  «  fie). 


Perturbation  State 


In  the  nondissipative  two-fluid  approximation  the  following  linearized  field  equations  describing 
the  perturbed  state  are  obtained  (s  =  e,  i). 

Equations  of  Momentum,  Mass,  and  Energy  Density  Convervation: 


—  +  v  ia  .  V 


vs  =  -  V  ps  + 


ns  es  [Eo 
nos  es[E+vsXBo+vasX  B  j 

vjns  =  -  (vs  .  V  nso  +  nos  V  •  VS) 

+  'is  •  V  j  PS  =  -  (Vg  .  V  PgQ  +  7g  Pg0  V  .  Vg) 

Maxwell's  Equations: 

1  X  B  =  ^  o  ^  (ng  eg  v  1S  +  nog  es  vg) 
s=e,  i 

VXE  = - 

at 


4  V  xs  X  Bo 


(69) 

(70) 

(71) 

(72) 

(73) 
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vm 


L/SON 


V  •  B  =  0 


V  •  E  =  * 


o'1  £  "s 


To  obtain  the  convective  solution  («~VqG)  of  Equations  (69)  through  (75),  any  equation  terms 
being  relatively  small  of  order  «  [(Equation  (43)]  are  neglected. 

Multiplication  of  Equation  (69)  by  (vO  gives,  after  neglecting  terms  ana  (v*  E^—**) 

'  7PS  +  nos  es  [e  +  vs  X  B0J  -  0  (76) 

In  the  same  approximation 


VXB=m0  £  noses?3 


is  obtained  from  Equation  (72). 

From  Equations  (77)  and  (76)  (s  *  e,  i)  results 

-w  - *  Bn  1  — •  -d  • 

-  7  (Pe  +  Pi)  x  B0  +  —  V  X  B - (B0  •  V  X  B)  B0  =  0  (78) 

Fo  ^o 

By  means  of  Equations  (76)  and  (78),  it  can  be  shown  from  Equation  (7  5)  that  the  plasma  behaves 
neutral  against  perturbations  with  wavelengths  X»d,  where  d  =  Jt Q  kTQ/noee^  is  the  Debye 
shielding  distance:30  ’ 

Z  ns  es  =  °  (79) 


Summation  of  Equations  (69)  (s  =  e,  i)  and  multiplication  by  V  X  gives,  because  of  Equations  (79) 
and  (72),  terms~t^  being  neglected: 


Z  nos  ms  ^,+vis'v 


~  1  X  1  *  _  1  _  S 

=  - (B0  •  V  )  V  X  B  + -  VX(B-VB0)+ - V  (B0  •  V  )  X  B 

F0  F  o  i - 1  F0 
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With  exception  of  0Bo  .  V  )  V  the  remaining  terms  in  Equation  (80)  are-*-*1.  In  apply¬ 

ing  the  €  -approximation  further,  Equation  (80)  predicts  that  the  convective  solution  is  associated 
with  perturbations  propagating  perpendicular  to  B0.  Now,  multiplying  Equation  (80)  by  7  X  gives 


(81) 


In  deriving  Equation  (81)  it  has  been  considered  that  k  X  B0,  and  v  .  0,  the  latter  following 

from  Equation  (70)  or  (71): 


[-|^-  +  V  ±B  •  vj  ps  +7s  •  V  Pso  =  -rs  Pso  V-  vg  ~  o  (82 

Eliminating  from  Equation  (81)  ve  by  means  of  Equation  (77),  then  VX§  by  means  of  Equation 
(78)  gives  (p  =  p.  +  fn 


and,  further,  after  multiplication  with  V  pD»  (7P0  -L  B0,  7  p0i-8  *  7  BQ)  and  by  consideration 
of  Equations  (64)  and  (66), 


Equation  (84)  becomes,  with  Equation  (82),  setting  there  s  =  i,  a  characteristics  equation  for 
the  ion  pressure  perturbation  p^ 


83 


From  Equation  (85)  one  derives  the  dispersion  equation: 

.2  (i  +n°ejM+1  2  .  r-  (v±i  - 

\  "oi  mi  /  V  noi  mt  x  j 

dT.  vxi)‘ 


‘oi  ‘“l 


(k.  v±er  -  o 


(86) 


Taking  into  account  the  relation  in  Equation  (67),  this  expression  reduces  to: 


According  to  Equation  (87),  the  particle  streams  are  unstable.  The  instability  is  caused  by  the 
inertia  forces  of  the  perturbed  electron  and  ion  streams,  which  are  in  interaction  with  each 
other  through  the  magnetic  field.  Both  instabilities  propagating  in  an  opposite  to  the  electron 

m 

drift  direction  are  excited.  The  effect  vanishes  for — —  —0. 

mi 


CONCLUSION 

Even  beyond  the  critical  Reynolds  number  above  which  the  ordinary  viscous  flow  instability 
sets  in,  magnetoactive  plasma  flows  can  become  unstable  due  to  convective  processes.  As  a 
result  of  these  instabilities,  the  directed  kinetic  energy  of  the  plasma  flow  is  partially  trans¬ 
formed  into  irregular  kinetic  and  potential  energy  of  nonequilibrium  fluctuations.  From  the 
large-scale  fluctuations,  this  kinetic  energy  is  transferred  by  nonlinear  interactions  to  the 
small-scale  fluctuations  and  by  the  latter  into  heat  (the  dissipation  mechanisms  in  the  plasma 
are  proportional  to  the  square  of  the  wave  number). 

With  respect  to  experimentation,  the  problem  becomes  that  of  determining  the  onset  and 
existence  of  turbulence  in  the  plasma  flow.  This  can  be  done,  for  example,  by  measuring  the 
intensity  variation  of  electromagnetic  waves  scattered  by  the  turbulence  elements.  Another 
method  which  can  be  proposed  is  to  add  strongly  radiating  test  atoms  to  the  plasma  at  one  side 
of  the  channel.  Assume  that  the  transverse  diffusion  time  of  the  test  atoms  is  large  with  re¬ 
spect  to  the  time  interval  within  which  the  test  atoms  are  carried  by  the  longitudinal  mass  flow 
through  the  channel.  Then,  in  laminar  plasma  flow  the  test  radiation  centers  cannot  be  ob¬ 
served  at  the  other  side  of  the  channel.  In  the  case  of  turbulent  plasma  flow,  however,  transverse 
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mass  motions  also  take  place.  In  the  latter  case,  part  of  the  test  radiation  centers  will  there¬ 
fore  reach  the  other  side  of  the  channel.  Their  concentration  can  then  be  determined  there  by 
measuring  the  intensity  of  the  characteristic  spectral  lines  of  the  test  atoms. 

In  the  state  of  developed  instability,  the  plasma  shows  anomalous  transport  properties.  It  is 
of  particular  interest  —  if  we  extend  the  results  to  partially  ionized  gases — how  the  coefficients 
of  recombination,  ionization,  and  the  ionization  potential  are  changed  by  the  nonequilibrium 
fluctuations.  The  investigation  of  the  latter  effects,  which  are  presumably  changing  the 
ionization  degree,  will  be  investigated  later. 
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